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ABSTRACT 


Whole  rock  Pb-Pb  and  Rb-Sr  isotopic  results  indicate 
an  age  of  deposition  of  less  than  2400-2470  m.y.  for  the 
sedimentary  precursors  of  the  Charlebois  Lake  complex. 

The  data  also  indicate  a  major  metamorphism  of  the 
precursors  during  the  Hudsonian  metamorphic  event  which, 
in  this  area,  ended  1750  m.y.  ago.  U-Th-Pb  determinations 
carried  out  on  oogenetic  uraninite  and  monazite-rich 
mineral  fractions  from  granodioritic  granofels  and  migma- 
tite  members  of  the  Charlebois  Lake  complex  establish 
open  system  chemical  behaviour  1750-1797  m.y.  ago.  Taken 
as  a  whole,  the  geochronological  evidence  suggests  the 
uranium  mineralization  is  a  result  of  the  local  remobili¬ 
zation  of  uranium-rich  sediments  during  the  Hudsonian 
metamorphic  event.  The  geochronological  conclusions 
agree  well  with  the  geological  conclusions  of  Morra  (1977) . 
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CHAPTER  I 


INTRODUCTION 

Because  of  sharply  rising  prices  of  uranium  over  the 
last  20  to  30  years,  low  grade  uranium  deposits  of  pegma¬ 
tite  type  have  become  increasingly  attractive  to  uranium 
exploration  companies. 

This  type  of  deposit  is  found  in  Precambrian  Shield 
areas  throughout  the  world,  most  commonly  occurring  in  the 
regionally  metamorphosed  terrains  that  have  been  exten¬ 
sively  granitized.  The  Canadian  Shield  contains  a  large 
number  of  these  sub-economic  deposits  and  although  no 
serious  exploitation  has  been  attempted  so  far,  these 
deposits  continue  to  be  studied  and  assessed  as  potential 
uranium  producers  in  the  future. 

One  of  the  better  known  U-pegmatite  deposits  in  the 
Churchill  province  is  in  the  Charlebois  Lake  area  in 
northeastern  Saskatchewan.  The  geology  of  the  area  has 
been  studied  in  detail  on  a  few  occasions  in  the  past  and 
therefore  it  is  considered  that  relatively  complete 
geological  documentation  for  the  area  is  available. 

The  purpose  of  this  study  was  to  investigate  the 
chronology  and  evolution  of  the  Charlebois  Lake  area  by 
examining  U-Th-Pb,  Pb-Pb  and  Rb-Sr  isotope  systems  from 
various  rock  units  of  the  Charlebois  Lake  complex.  It  was 
hoped  that  these  isotopic  systems  would  record  events  that 
affected  the  area  and  would  provide  insight  into  the 
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origin  of  the  Charlebois  Lake  complex  and  associated  U- 
mineralization . 

The  samples  were  collected  so  that  they  represent  the 
major  units  of  the  complex:  granitic  to  tonalitic  gneiss, 
mineralized  granodioritic  granofels,  biotite  and  hornblende 
gneiss,  amphibolite  and  late  felsic  dyke. 


CHAPTER  II 


REGIONAL  AND  LOCAL  GEOLOGY 

The  Charlebois  Lake  area  is  located  in  northeastern 
Saskatchewan,  about  160  km  west  of  the  Manitoba  border 
and  50  km  south  of  the  Northwest  Territories  border. 
Approximate  coordinates  are  59°50'N  and  104°40'W 
(Figure  1) .  Although  the  geological  surveys  of  northern 
Saskatchewan  were  initiated  at  the  beginning  of  this 
century  and  intensified  after  the  Second  World  War,  with 
the  increasing  need  for  uranium  most  of  the  efforts  were 
concentrated  primarily  north  of  Lake  Athabasca  and  west 
of  the  Charlebois  Lake  area. 

The  first  serious  reconnaissance  work  in  the  area 
was  done  by  Mawdsley  (1950)  who  described  a  number  of 
radioactive  occurrences.  This  work  was  followed  by  more 
detailed  field  work,  mapping,  petrological  and  petrogra¬ 
phic  studies  by  Mawdsley  (1952),  Cumming  (1952)  and 
Kirkland  (1952) .  From  that  time  on,  for  a  period  of 
about  twenty  years,  not  much  activity  has  been  reported 
in  this  area. 

In  the  meantime  regional  work  on  the  Saskatchewan 
Shield  continued.  Geophysical  surveys  of  northern 
Saskatchewan  were  carried  out.  The  gravity  and  aeromag- 
netic  results  suggested  that  the  Saskatchewan  crystalline 
basement  was  composed  of  several  distinct  tectonic  units 
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3.  Needle  Tails  Shear  zone.  4.  Parker  L  Shear  /one 
5.  Tabbernor  Fault 

/  Otfier  I  it  host  rue  t  era  I  domain  boundaries 

Figure  1.  Lithostructural  subdivision  of  Northern  Saskatchewan 
(from  Lewrey  et  al.  1978) . 


. 


'V-  -V  '  -  ,t&.  I  r  i  \  v 


5 


which  are  traceable  to  the  Cordillerian  fold  belt.  As 
such  they  were  perceived  as  deep  crustal  entities.  Further 
geological  studies  of  the  Saskatchewan  Precambrian  base¬ 
ment  were  of  broader  scope  and  included  structural,  petro¬ 
logical,  geochemical  and  also  a  number  of  geochronological 
investigations . 

Four  lithos tructural  units  presently  recognized  in 
Saskatchewan  were  discussed  and  described  in  a  number  of 
important  papers,  particularly  those  by  Lewrey  et  al. 
(1978),  Burwash  (1978,1979),  Sibbald  et  al .  (1976), 

Burwash  et  al.  (1962,1969,1970),  and  Foster  et  al.  (1970). 

A  few  of  these  publications  also  dealt  with  main  events 
leading  to  the  final  stabilization  of  this  part  of  the 
Churchill  structural  province.  As  a  result  it  is  now 
certain  that  much  of  the  Archean  craton  in  Saskatchewan 
was  reactivated  during  the  Hudsonian  orogeny  which  lasted 
from  1900  m.y.  ago  to  1700  m.y.  ago  (Burwash  et  al.  1969). 
Reactivation  consisting  of  extensive  granit i zation ,  K- 
metasomatism  and  polymetamorphism  affected  much  of  the 
Saskatchewan  Shield. 

The  four  lithostructural  units  of  the  Saskatchewan 
Shield  are  shown  in  Figure  1.  Significant  features  of 
these  northeasterly  trending  units  are  described  in  detail 
in  the  paper  by  Lewrey  et  al.  (1978).  A  brief  summary  of 
two  of  these  units  will  be  given  here;  since  the 
Charlebois  Lake  area  is  part  of  one  unit,  it  is  geogra- 
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phically  close  to  the  other  and  together  they  contain 
most  of  the  radioactive  anomalies  discovered  so  far  in 
Saskatchewan . 

The  first  is  the  Western  Craton  which  occupies  the 
northeast  corner  of  Alberta  and  western  Saskatchewan. 

It  is  comprised  of  an  Archean  basement  complex  of  grani¬ 
toids,  migmatites  and  other  high  grade  metamorphic  rocks. 
Reactivation  of  the  Craton  during  the  Hudsonian  orogeny 
resulted  in  a  series  of  younger  mobile  belts  which 
surround  the  older  stable  blocks  (Beck  1969).  Hudsonian 
overprinting  occurred  under  conditions  of  lower  amphi¬ 
bolite  and  greenschist  facies  metamorphism.  Here,  depend¬ 
ing  on  the  degree  of  reactivation,  a  range  of  ages  from 
2500  m.y.  to  1750  m.y.  can  be  observed  (Baadsgaard  and 
Godfrey  1970,  Sassano  et  al.  1972,  Koeppel  1968). 

In  this  area  a  large  number  of  uraniferous  pegmatite 
deposits  can  be  found.  Most  of  them  are  pegmatite  deposits 
of  group  A  (Trembley  1978) .  One  of  the  more  prominent  is 
the  Grease  River  occurrence  which  is  located  only  about 
30  km  west  of  the  Charlebois  Lake  area.  Here  uranium 
occurs  in  massive  to  faintly  foliated  granite  and  in  the 
pegmatites  that  form  part  of  the  Archean  basement  gneiss 
complex.  The  uraniferous  pegmatites  are  of  various 
dimensions  and  are  erratically  distributed.  Since  the 
pegmatites  occur  only  in  the  basement  granites  and 
granitic  gneisses  to  which  they  are  closely  related  in 
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composition,  they  are  regarded  as  the  anatectic  fractions 
of  the  basement  granitic  rocks.  They  are  generally  coarse¬ 
grained  and  not  related  to  a  single  large  intrusive  granite 
body . 

The  central  part  of  the  Saskatchewan  Shield  is 
occupied  by  the  Cree  Lake  zone.  This  zone  is  flanked  to 
the  west  by  the  Virgin  River  domain  and  to  the  east  by  the 
Wollaston  domain.  The  zone  is  comprised  of  mainly  Archean 
granitoids,  migmatites  and  high  grade  metamorphic  rocks. 

The  whole  area  is  characterized  by  foliated,  folded  and 
boudinaged  small  scale  granitoid  segregations  which  might 
indicate  extensive  "in  situ"  anatexis.  It  is  suggested 
that  pressure-temperature  conditions  throughout  most  of 
the  Cree  Lake  zone  during  the  metamorphism  were  at  least 
as  high  as  those  corresponding  to  the  upper  amphibolite 
facies . 

Lewrey  et  al.  (1978)  state  that  in  the  Cree  Lake  zone 
"it  is  probable  that  the  high  grade  conditions  achieved 
early  in  the  def ormational  history  were  maintained  at 
essentially  the  same  level  throughout  and  beyond  the  close 
of  the  major  deforma tional  event.  Accompanying  anatexis 
was  effectively  'in  situ'  as  there  is  little  evidence  of 
large  scale  intrusion". 

Some  reference  to  this  statement  will  be  made  again 
in  Chapter  III  when  discussing  the  ages  of  uraninites  and 
monazites.  Metamorphic  grade  of  the  zone  falls  slightly 
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in  the  cooler,  more  brittle,  marginal  areas  of  Virgin 
River  and  Wollaston  domain  (Figure  2) . 

Geochronological  studies,  although  limited  mainly  to 
the  Wollaston  domain  and  areas  south  of  Athabasca  basin, 
yield  a  range  of  ages  from  2700  m.y.  to  1570  m.y.  (Cumming 
and  Scott  1975,  Weber  et  al .  1975,  Money  et  al .  1970)  . 

Uraniferous  pegmatite  deposits  are  numerous  in  the 
Cree  Lake  zone.  Both  Group  A  and  Group  B  deposits  can  be 
found  (Trembley  1978).  Group  B  deposits  (Charlebois  Lake 
occurrence  is  a  typical  representative  of  this  deposit) 
are  characterized  by  the  concordance  of  the  pegmatite 
layer  with  stratigraphically  lower  granitic  gneisses  and 
stratigraphically  higher  mafic-rich  gneisses,  gradational 
contacts  between  the  pegmatite  and  other  rock  units  and 
much  more  variable  grain  size  and  composition  than 
exhibited  by  Group  A. 

In  1974,  Fosago  Exploration  Ltd.  carried  out  an 
extensive  exploration  program  in  the  Char lebois-Higginson 
Lake  area,  consisting  of  detailed  geological  mapping  and 
radioactive  surveys  in  order  to  study  the  relative  strati¬ 
graphic  position  of  the  radioactive  occurrences.  In  1977 
F.  Morra  completed  an  M.Sc.  thesis  in  which  he  described 
the  geology  and  uranium  deposits  of  the  Charlebois  Lake 
area.  This  thesis  (along  with  the  work  of  Cumming  (1952), 
Kirkland  (1952)  and  Mawdsley  (1952))  was  used  as  a  major 
source  of  geological  reference  in  the  present  study. 
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Figure  2.  Metamorphic  map  of  the  Saskatchewan  Shield  (from  Lewrey  et  al.  1978)  . 


10 


Most  of  the  Charlebois  Lake  area  which  is  situated 
in  the  northwestern  part  of  the  Cree  Lake  zone  and  which 
appears  to  be  the  northern  extension  of  the  Virgin  River 
domain,  is  underlain  by  massive  and  foliated  granitic  to 
tonalitic  gneiss.  Other  lesser  rock  units  are  migmatite, 
granodiorit ic  granofels  (or  younger  intrusive  pegmatite, 
as  called  by  Cumming,  1952),  calc-silicate  rocks,  horn¬ 
blende  gneiss,  amphibolite,  biotite  gneiss  and  quartzite. 
Figure  3  (a  map  adopted  from  Morra  (1977))  gives  the 
relative  distribution  of  these  rock  units. 

All  of  the  rock  units  in  the  area  are  gently  to 
strongly  folded  with  few  poorly-exposed  faults.  The  most 
prominent  lineament  is  represented  by  the  Higginson  Lake 
fault  which  is  believed  to  be  an  offset  extension  of  the 
regionally  important  Clut  Lake  fault.  North  of  this  fault 
there  are  possibly  a  few  more  parallel  faults:  one 
through  Dramnitzke  Bay,  one  immediately  north  of  it,  one 
from  Heisler  Lake  to  Sprekley  Lake,  and  others.  One  of 
the  rare  occurrences  of  the  mylonite  zone  has  been  observed 
by  Cumming  (1952)  on  the  granitic  gneiss  on  the  west  shore 
of  Sprekley  Lake.  This  area  is  presumably  an  extension  of 
the  east-west  trending  fault  that  originates  at  Heisler 
Lake . 

The  largest  number  of  folds  follows  the  general 
regional  and  local  fault  trends,  i.e.  axes  of  the  folds 
have  a  northeast-southwest  direction.  The  cores  of  the 
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Figure  3.  Geological  map  of  the  Charlebois  Lake  Area  (from  Morra,  1977) 


anticlines  are  occupied  by  the  granitic  to  tonalitic 
gneisses  while  synforms  are  generally  occupied  by  other 
rock  units,  such  as  biotite,  hornblende  gneisses,  calc- 
silicates  and  quartzites. 

The  granitic  to  tonalitic  gneiss  is  best  represented 
south  of  Dramnitzke  Bay  where  it  forms  the  core  of  the 
overturned  Pegasus  Lake  anticline.  The  anticline  is 
flanked  on  both  sides  by  a  metasedimentary  sequence  which 
is  continuous  and  regular. 

In  an  unpublished  report,  Sassano  (1974)  attempted  to 
correlate  the  granitic  to  tonalitic  gneiss  of  this  area 
with  the  Donaldson  Lake  gneiss  (of  Middle  Tazin  age)  from 
the  Beaverlodge  area  (Sassano  1972).  He  also  suggested 
that  an  unconformity  separates  the  granitic  and  tonalitic 
gneiss  from  the  rest  of  the  metasedimentary  sequence, 
which  he  tentatively  correlated  to  the  Fay  Mine  sequence 
(of  Upper  Tazin  age) .  Cumming  (1952)  suggested  that  the 
coarse  to  very  coarse  grained,  mineralized  rock  which 
outcrops  around  Pegasus  Lake  granitic  gneiss  is  in  fact 
younger  intrusive  pegmatite. 

Based  on  field  observations,  stratigraphic  and 
structural  characteristics  such  as  concordance  of  the 
contacts,  regularity  of  the  sequence  of  the  rock  units  in 
the  antiforms  and  synforms,  petrographic  and  other  labora¬ 
tory  studies,  Morra  (1977)  concluded  that  "young  intrusive 
pegmatite"  (which  he  calls  granodior itic  granofels)  is  a 


part  of  the  metasedimentary  sequence.  He  also  concluded 
that  the  granitic  to  tonalitic  gneiss  is  the  oldest,  basal 
and  integral  part  of  the  Charlebois  Lake  metasedimentary 
complex . 

A  list  of  stratigraphic  units,  as  summarized  by  Morra 
(1977),  is  given  in  Table  I. 

Since  the  lower  contact  of  granitic  and  tonalitic 
gneiss  is  not  visible  in  the  area,  it  is  believed  that  the 
Charlebois  Lake  complex  could  have  a  thickness  of  up  to 
3000  meters  and  that  the  crystalline  basement  rocks  are 
not  exposed. 

The  rocks  of  the  area  have  undergone  regional  meta¬ 
morphism  in  the  intermediate  amphibolite  facies.  The 
presence  of  sillimanite  and  garnet  associated  with  biotite, 
microcline  and  oliogoclase,  along  with  other  metamorphic 
assemblages,  is  widespread  (Cumming  1952,  Morra  1977). 

This  paragenesis  was  described  by  Winkler  (1967)  as  the 
sillimanite-almandine-orthoclase  subfacies  of  the  amphi¬ 
bolite  facies.  Mineral  assemblages  belonging  to  the 
upper  limit  of  amphibolite  facies,  although  not  very 
common,  were  also  found  by  Morra  (1977).  This  seems  to  be 
in  good  agreement  with  regional  metamorphic  patterns. 
Lewrey  et  al .  (1978)  suggested  tha.t  upper  amphibolite 

conditions  have  persisted  throughout  the  Cree  Lake  zone 
with  a  very  gentle  decrease  in  metamorphic  grade  towards 
the  eastern  and  western  margins  of  the  zone.  The  grade 
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TABLE  I 


TABLE  OF  STRATIGRAPHIC  UNITS  (FROM  MORRA  1977) 
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Athabasca  Formation  (south  of 
the  mapped  area) . 
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Martin  Formation  (west  and 
possibly  south  of  the  mapped 
area) . 

Major  Unconformity 


Pink  pegmatite  and  irregular 
bodies 


Intrusive  contact 


Quartzite 
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Biotite  Gneiss 
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Hornblende  Gneiss  and 
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Aphebian 
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Calc-silicates 

VI 
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O 

Migmatite 
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w 

Granodiorit ic  Granofels 

pp 

< 

na 

u 

Granitic  and  Tonalitic 

_ 

Gneiss 
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decreases  very  rapidly  within  a  few  kilometers  of  the 
margins.  In  the  Virgin  River  Shear  Zone  a  decrease  from 
granulite  facies  to  lower  amphibolite  facies  takes  place 
within  about  13  kilometers.  Although  the  situation  in 
the  northern  part  of  the  Cree  Lake  zone  is  probably  not 
the  same  as  in  the  areas  south  of  the  Athabasca  Basin, 
there  are  some  indications  that  the  metamorphic  grade 
between  Charlebois  Lake  and  the  Black  Lake  fault  decreases. 
Johnston  (1963,1964)  has  observed  a  large  volume  of  schists 
in  the  Lytle  Lake  area.  The  Grease  River  area  (Costaschuk 
1979)  which  is  about  20  km  west  of  the  Black  Lake  fault, 
has  undergone  metamorphism  in  lower  amphibolite  to  green- 
schist  facies.  The  Grease  River  area  is  a  part  of  the 
Western  Craton  and  not  of  the  Cree  Lake  zone  so  this 
comparison  might  not  be  appropriate.  Unfortunately  there 
is  not  enough  published  evidence  to  substantiate  or  dis¬ 
prove  similarities  between  the  marginal  areas  of  the  Cree 
Lake  zone  south  and  north  of  the  Athabasca  sandstone  basin. 

Based  on  the  pronounced  structural  concordance  of  the 
members  of  the  Charlebois  Lake  complex,  gradational 
contacts,  strong  foliation  and  banding,  continuity  of  the 
units  and  the  presence  of  conformably  intermixed  calc- 
silicates  with  all  of  the  other  upits,  Morra  (1977) 
suggested  that  pre-existing  rocks  were  sediments  deposited 
in  dry  and  shallow  marine  environments  which  were  subjected 
to  sudden  and  cyclic  changes  during  the  deposition.  The 
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original  sedimentary  sequence  consisted  of  arkosic  are- 
nites ,  calcareous  pelites,  sandstones  and  tuffaceous  and 
pyroclastic  material,  pelites  and  feldspathic  arenites. 

The  present  rocks  are  mono-metamorphic  derivatives  of 
these  sediments.  The  metamorphic  potash  migration,  which 
has  not  been  on  a  very  large  scale,  might  have  been  partly 
responsible  for  migmatite  formation. 

The  late  pink  felsic  dykes  and  irregular  pegmatitic 
bodies  outcrop  at  various  locations,  crosscutting  almost 
all  of  the  other  rock  units.  According  to  Cumming  (1952) 
and  particularly  Morra  (1977),  the  contacts  between  these 
late  rocks  and  the  rest  of  the  units  of  the  Charlebois  Lake 
complex  appear  to  be  the  only  unconf ormable ,  intrusive 
contacts  in  the  area. 

No  detailed  geochronological  studies  have  been 
carried  out  so  far  in  the  Charlebois  Lake  or  immediately 
surrounding  areas.  The  only  published  isotopic  data  from 
the  Charlebois  Lake  area  are  by  Cumming  et  al.  (1955), 
obtained  on  uranium  minerals.  Recalculation  of  original 
data  yields  a  207Pb/206Pb  age  of  1765  m.y.  and  1785  m.y. 
for  two  uranium  minerals. 

Detailed  mineralogical  and  petrological  descriptions 
of  all  of  the  rock  units  are  given  by  Cumming  (1952), 
Kirkland  (1952)  and  Morra  (1977) .  What  follows  is  only 
a  brief  summary  of  the  main  features  of  each  of  the  rock 


units . 
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Granitic  and  tonalitic  gneiss 

This  rock  unit  represents  more  than  50%  of  the 
exposed  rocks  in  the  area.  Initially  it  was  described 
and  mapped  as  "late  intrusive  formation"  or  orthogneiss 
(Mawdsley  1952)  and  subsequently  termed:  granite,  pegma¬ 
tite  granite,  granite  gneiss,  massive  granite  etc.  It  is 
usually  exposed  in  anticline  areas  of  which  a  very  good 
example  is  the  overturned  Pegasus  Lake  anticline. 

The  rock  grades  from  massive,  fine-grained,  gray 
gneiss  in  the  center  of  the  anticline  to  slightly  foliated 
in  the  outer  margins,  to  strongly  foliated  (when  the 
biotite  content  is  larger  than  8%)  near  the  contacts 
with  the  other  units. 

The  gneiss  is  pinkish  in  colour,  in  the  more  foliated 
type,  and  generally  K-feldspar  content  is  higher.  Grain 
size  varies  from  0.5  mm  to  1  mm.  Concordant  amphibolite 
bands  of  up  to  2  meters  thick  are  present  throughout  the 
gneiss  but  are  especially  prevalent  near  the  contact  with 
the  other  units. 

The  mineralogy  is  rather  monotonous.  Quartz,  plagio- 
clase,  K-feldspar  and  biotite  are  the  main  constituents. 

The  relative  proportions  of  these  minerals  changes  depend¬ 
ing  on  the  degree  of  foliation. 

Granodior itic  granofels 

"Pegmatite  granite"  (Mawdsley  1950) ,  "young  intrusive 
fine  pegmatite"  (Cumming  1952)  and  "granodioritic  granofels" 
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(Morra  1977)  are  common  terms  used  to  describe  this  rock 
unit . 

It  is  well  represented  along  the  southern  shore  of 
Dramnitzke  Bay,  having  a  thickness  of  30  m  to  90  m.  It 
borders  the  granitic  gneiss  as  a  continuous  lenticular 
body  which  is  parallel  to  the  regional  foliation.  Horizons 
of  granodioritic  granofels  often  alternate  with  the  hori¬ 
zons  of  migmatite.  The  contacts  are  gradational  and 
conformable . 

Morra  (1977)  observes:  "Granodioritic  granofels  does 
not  demonstrate  a  crosscutting  (or  intrusive)  relationship 
with  the  adjacent  units  anywhere  in  the  area.  Its  paralle¬ 
lism  with  the  stratigraphically  lower  granitic  to  tonalitic 
gneiss  and  higher  calc-silicate  is  constant  and  continuous". 
He  suggests  that  granofels  is  the  most  appropriate  term 
for  this  rock  since  it  is  compatible  with  the  textural 
and  mineralogical  characteristics  and  it  does  not  imply 
that  the  rock  was  either  igneous  or  sedimentary  before  the 
metamorphism,  while  "pegmatite"  carries  a  genetic  implica¬ 
tion  . 

Mineralogically ,  the  rock  is  indistinguishable  from 
a  granitic  gneiss.  Grain  size  is  the  main  distinguishing 
feature  ( 3  mm  to  30  mm  compared  to  0.5  mm  to  1  mm  for 
granitic  gneiss) .  Granodioritic  granofels  show  very  little 
foliation  but  the  grain  size  becomes  finer  and  biotite 
assumes  subparallel  orientation  as  one  approaches  the 
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contact  with  the  granitic  gneiss.  The  colour  is  light 
gray  to  white,  and  occasionally  dark  gray  when  biotite  is 
very  abundant.  Besides  quart z,  plagioclase ,  K-feldspar 
and  biotite,  the  rock  commonly  contains  variable  amounts 
of  uraninite  and  secondary  uranium  minerals,  monazite, 
zircon,  sphene  and  molybdenite. 

Migmatite 

The  migmatites  are  composed  of  biotite  rich  portions 
and  separated  lenses  of  granitic,  granodior itic  and  tona- 
litic  composition.  The  formation  of  this  rock  unit  could 
be  due  to  relatively  limited  potash  migration  during  the 
metamorphic  event  or  it  could  be  due  to  a  rather  gradual 
and  cyclic  change  in  the  depositional  environment  of  the 
original  sediments.  Its  mineralogical  composition  is 
intermediate  between  two  rock  types:  quartz-f eldspathic 
bands  (composition  equivalent  to  granitic  gneiss)  and 
biotite  rich  bands  (equivalent  to  the  biotite  gneiss). 

It  is  suggested  (Morra  1977)  that  in  terms  of  original 
deposition,  migmatites  correspond  to  the  transitional 
period  of  arkosic  and  pelitic  sedimentation,  repeated  in 
short,  distinct  cycles,  which  gave  rise  to  conformable  and 
regularly  interbedded  felsic  (arkosic)  and  mafic  (pelitic) 
bands . 

The  felsic  intercalating  bands  are  mainly  composed 
of  quartz,  plagioclase  and  K-feldspar,  in  various  propor¬ 
tions.  The  mafic  bands  are  biotite-and/or  hornblende-rich, 
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are  well  foliated  and  are  usually  finer  grained  than 
felsic  bands  (approximately  0.5  mm  and  up  to  4  mm,  res¬ 
pectively)  .  Garnet  rich  bands  are  present  and  serve  as 
marker-horizons  for  this  unit.  Uraninites  and  molyb¬ 
denites  are  common,  as  well  as  other  accessory  minerals 
such  as  apatite,  sphene,  sillimanite,  and  occasionally 
pyrite  and  monazite. 

Calc-Silicates 

This  unit  includes  amphibole-diopside  gneiss,  talc- 
serpentine-carbonate  marble  and  ser icite-phlogopite/ 
biotite-diopside  gneiss.  The  horizon  of  this  unit, 
which  stratigraphically  occurs  between  granodior it ic 
granofels  and  hornblende  gneiss,  is  present  on  the  shore 
of  Charlebois  Lake,  Sprekley  Lake  and  Dramnitzke  Bay. 

These  "softer"  gneisses  have  been  gouged  away  by  glaciers 
and  are  generally  extensively  weathered.  This  rock  unit 
could  be  up  to  400  m  thick,  appears  as  a  continuous  band 
for  several  miles  and  follows  the  general  foliation 
trend  of  the  area. 

Grain  size  is  fine  to  medium,  it  generally  exhibits 
a  massive  habit  and  is  dirty-white  to  gray-green  in  colour. 

Hornblende  gneiss  and  amphibolite 

This  unit  is  best  represented  north  of  Dramnitzke 
Bay  where  it  is  closely  associated  with  biotite  gneiss  in 
the  northeast-southwest  trending  belt.  It  contains  an 
increasing  amount  of  biotite  towards  the  contact  with 
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biotite  gneiss,  suggesting  a  gradational  character  between 
these  two  units. 

The  colour  of  the  rock  depends  on  the  amount  of  horn¬ 
blende  and  varies  from  medium  gray  to  dark  gray.  The  grain 
size  is  usually  fine  to  very  fine  (less  than  0.5  mm)  and 
habit  is  massive  if  the  amount  of  hornblende  is  more  than 
50%  and  little  quartz  is  present.  K-feldspar  is  generally 
absent  or  present  in  small  amounts. 

Biotite -gneiss 

This  gneiss  is  also  well  represented  in  the  area 
immediately  north  and  west  of  Dramnitzke  Bay  in  the 
northeast-southwest  striking  belt,  and  is  intimately 
associated  with  the  hornblende  gneiss.  The  percentage 
of  hornblende  increases  towards  the  contact  with  horn¬ 
blende  gneiss.  Morra  (1977)  states:  "This  unit  is 

difficult  to  map  because  of  the  presence  of  numerous 
granodioritic  granofels  and  hornblende  gneiss  intercala¬ 
tions.  Therefore  most  of  the  area  mapped  as  biotite 
gneiss  may  also  contain  hornblende  gneiss  as  well  as 
3-6  meters  thick  granodioritic  granofels  bands".  The 
rock  is  strongly  foliated  due  to  the  alignment  of  biotite 
grains.  It  is  medium  gray  in  colour  and  grain  size  is 
fine  to  medium,  averaging  0.5  mm.  The  main  constituent 
minerals  are  quartz,  plagioclase,  microcline  and  biotite. 
Accessories  are  metamict  zircon,  apatite  and  muscovite, 
and  occasionally  sillimanite  and  sphene. 


Quartzite 


Strat igraphically  this  unit  should  be  the  youngest 
member  of  the  Charlebois  Lake  complex.  It  is  found  as 
small  occurrences,  beds  or  lenses  within  the  area, 
generally  following  conformable  biotite  gneiss.  The 
quartzite  is  not  "pure".  Quartz  content  averages  75%, 
the  remaining  portion  being  feldspar  with  a  minor  amount 
of  biotite. 

Late  felsic  dykes  and  irregular  bodies 

All  of  the  units  of  the  Charlebois  Lake  complex  are 
crosscut  by  these  late  felsic  dykes  at  different  loca¬ 
lities.  The  dykes  are  generally  3-6  meters  thick  with 
the  grain  size  ranging  from  medium  to  coarse,  while 
irregularly  shaped  pegmatite  bodies  occupy  substantially 
larger  areas  and  the  grain  size  is  coarse  to  very  coarse 
(individual  grains  occasionally  reaching  10  cm  in 
diameter) .  The  colour  is  salmon  pink  due  to  the  presence 
of  a  large  amount  of  microcline.  The  contact  with  the 
other  units  is  sharp  and  postcrystalline  deformation  is 
evident  in  many  areas.  Major  mineral  constituents  are 
microcline,  quartz,  some  plagioclase  and  very  small 


amounts  of  biotite. 


CHAPTER  III 


RESULTS  AND  DISCUSSION 

A  suite,  consisting  of  over  forty  rock  samples  was 
collected  during  two  field  trips.  First,  during  the 
summer  of  1976  by  Dr.  G.L.  Cumming  of  the  Physics  Depart¬ 
ment  and  second,  in  the  fall  of  1977,  by  Dr.  Cumming  and 
the  author. 

Extensive  weathering  of  some  of  the  rock  units, 
transportation  and  other  difficulties  did  not  permit 
collection  of  large  specimens  (by  volume) ,  nor  did  they 
allow  sampling  of  all  the  members  of  the  Charlebois  Lake 
complex.  The  size  of  the  samples  ranged  from  approxi¬ 
mately  one  pound,  for  the  specimens  collected  along  the 
southern  shore  and  north  of  Dramnitzke  Bay,  to  a  few 
pounds  for  the  specimens  from  the  center  of  the  Pegasus 
Lake  anticline. 

Rock  units  sampled  include  granitic  to  tonalitic 
gneiss,  granodioritic  granofels,  migmatite,  biotite  gneiss, 
hornblende  gneiss  and  late  felsic  dyke.  Sample  locations 
and  identification  numbers  are  shown  in  Figure  4. 

For  convenience  and  ease  of  interpretation,  all  of 
the  samples  are  divided  and  considered  as  three  groups: 
the  "Pegasus  Lake  Area"  group  which  includes  samples  1  to 
20,  is  represented  by  specimens  which  range  from  gray 
massive  granitic  to  tonalitic  gneisses,  to  more  foliated 
gray  to  pinkish  granitic  gneisses  and  three  specimens  of 
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Figure  4.  Sample  location  map. 
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coarse  grained,  mineralized  granodioritic-granof els . 

The  area  north  of  Dramnitzke  Bay  where  samples  21  to 
34  were  collected,  was  traditionally  the  field  camp  area 
of  Mawdsley  (1950),  Cumming  (1952),  Kirkland  (1952)  and 
Fosago  Exploration  (Morra  1977).  This  area  will  be 
referred  to  as  the  "Old  Camp  Area" .  The  area  is  comprised 
of  migmatites  along  the  shore  of  Dramnitzke  Bay,  which 
northward  grade  to  biotite  gneiss,  hornblende  gneiss, 
culminating  with  amphibolites,  some  three  to  five  hundred 
meters  from  the  shore  of  Dramnitzke  Bay.  A  not  very  well 
defined  fault,  following  the  regional  and  local  northeast- 
southwest  trend,  appears  to  be  present  in  the  hornblende 
gneiss-amphibolite  area.  Northwest  of  this  apparent  fault 
the  rocks  again  grade  from  gneiss  with  high  mafic  content 
to  almost  gray  granitic  gneiss. 

The  third  group  of  samples  (35  to  40)  consists  of 
specimens  collected  from  the  fifteen  foot  wide  late  pink 
felsic  dyke  which  crosscuts  the  other  rock  units  in  the 
"Old  Camp  Area". 

The  sample  size  and  the  fact  that  not  all  the  members 
of  the  Charlebois  Lake  complex  have  been  represented,  have 
put  definite  limits  on  the  depth  of  this  study,  as  will 
be  seen  later  in  this  chapter. 

Most  of  the  calculations,  including  isotopic  ratios, 
concentrations,  straight  line  fitting  (Cumming  et  al.  1972), 
and  the  drawing  of  diagrams  were  done  using  computer  pro¬ 
grams  written  by  Dr.  G.L.  Cumming.  Ages  and  error  limits 
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from  the  intersections  of  lines  with  the  concordia  or 
growth  curve  were  calculated  using  computer  programs 
written  by  the  author. 

Whenever  a  straight  line  was  fitted  to  a  set  of 
points,  the  errors  used  in  the  calculations  were  at  one 
sigma  level.  Assigned  correlation  coefficients  were: 

0.95  for  lead/uranium  data,  0.75  for  lead/lead  data  and 
0.0  for  rubidium/strontium  data.  Therefore,  throughout 
this  study,  all  the  slopes  are  given  at  one  sigma  error 
level.  All  the  age  values  obtained  from  the  slope  of  the 
line  or  from  the  intercepts  of  lines  with  the  concordia 
or  growth  curve,  were  calculated  at  two  sigma  error  level 
including  the  error  resulting  from  the  imperfect  fit  of 
the  data,  as  given  by  MSWD  (Mean  Squared  Weighted  Deviates) 
factor.  The  MSWD  factor  is  a  measure  of  how  well  the  data 
points  fit  the  line  within  the  assigned  analytical  error 
limits.  In  the  case  of  a  perfect  fit,  MSWD  will  have  a 
value  of  1.  Values  higher  than  1  suggest  that  the  data 
points  scatter  due  to  "geological"  reasons  (providing  the 
analytical  errors  have  not  been  underestimated) .  Values 
smaller  than  1  indicate  that  the  data  points  probably  fit 
the  line  but  that  the  analytical  error  limits  have  been 
overestimated.  A  very  brief  treatment  of  the  geochrono- 
logical  theory  applicable  to  the  isotopic  systems  studied 
in  this  thesis,  along  with  a  set  of  important  general 
references,  are  given  in  Appendix  I. 
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U-Th-Pb  Results 

Five  rock  specimens  were  available  for  this  work. 

They  included  three  samples  of  granodioritic  granofels 
(two  very  coarse  grained  and  one  medium  grained)  from  the 
mineralized  area  on  the  southern  shore  and  two  samples  of 
migmatite  with  high  whole  rock  lead-lead  ratios  from  the 
north  shore  of  Dramnitzke  Bay. 

The  rock  samples  were  crushed  and  pulverized.  Pb-Pb 
and  Rb-Sr  analyses  were  performed  on  the  whole  rock  powder, 
after  which  the  heavy  minerals  were  separated.  Small 
relatively  clean  uraninite  grains  and  monazite-r ich  frac¬ 
tions  were  picked  out  under  the  microscope  from  three 
specimens  (18,  19  and  22),  "uraninite"  from  specimen  20 
was  extracted  by  leaching  the  heavy  mineral  separate  with 
4N  HNO^  and  after  that  the  monazite-rich  phase  was  con¬ 
centrated.  Specimen  23  did  not  contain  any  visible 
uraninite  so  that  only  a  monazite-rich  fraction  was 
prepared.  One  specimen  of  coarse  grained,  biotite  rich, 
radioactive  pegmatite  from  Sickle  Lake  (southwest  of 
Pegasus  Lake  antiform) ,  collected  in  1951  and  analysed  at 
the  University  of  Toronto  (Cumming  1955)  was  also  avail¬ 
able.  A  polished  section  was  prepared  and  very  small 
grains  of  what  was  assumed  to  be  uranium  mineral  were 
separated  under  the  microscope  by  Dr.  R.D.  Morton  of 
this  department. 

The  samples  were  analysed  for  lead,  uranium  and  some 
for  thorium  isotopes.  Detailed  analytical  procedures  are 
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given  in  Appendix  II.  The  analytical  results  are 
presented  in  Tables  II,  III  and  IV.  Two  sigma  errors 
of  1%  were  assigned  to  the  207Pb/235U,  206Pb/238U  and 

208Pb/232Th  ratios,  as  suggested  in  Appendix  II.  Lead- 
uranium  ratios  are  plotted  on  the  standard  concordia 
diagram  (Figure  5).  The  size  of  the  elipses  represents 
the  error  limits  at  three  sigma  level  and  the  subscripts 
"U"  or  "M"  attached  to  the  sample  identification  number 
designates  a  uraninite  or  monazite-rich  phase. 

It  should  be  pointed  out  that  the  uranium,  thorium 
and  lead  concentration  results  are  to  be  viewed  as 
relative  rather  than  absolute  with  respect  to  the  analysed 
specimen.  This  is  particularly  so  where  the  total  avail¬ 
able  sample  is  very  small,  the  phases  are  impure  or  where 
leaching  was  involved. 

The  data  points  in  the  concordia  diagram  appear  to 
form  two  distinctly  different  trends.  One  uraninite 
sample  (22U)  is  concordant  and  plots  on  the  concordia 
curve.  Two  other  uraninite  separates  exhibit  normal 
discordance  (18U  and  19U) ,  i.e.,  they  plot  in  the  lead 
loss  field.  The  leached  uraninite  (20U)  exhibits  slightly 
reverse  discordance  (Baadsgaard  1965) .  Two  monazites 
plot  in  the  lead  loss  field  (18M  and  19M)  and  the  other 
three  show  reverse  discordance  (relatively  common 
behaviour  for  monazites) .  The  Sickle  Lake  specimen 
contains  about  30%  common  lead  and  shows  very  high  reverse 
discordance.  The  high  content  of  common  lead  makes  it 
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very  susceptible  to  error  due  to  the  uncertain  isotopic 
composition  of  common  lead,  so  that  any  interpretation 
of  U-Pb  data  becomes  highly  unreliable.  For  that  reason 
this  sample  has  not  been  plotted  on  the  concordia  diagram. 

Three  uraninite  data  points  (excluding  leached 
uraninite  20U)  form  one  group  and  the  monazite-rich  phases 
form  the  other.  Leached  uraninite  (20U)  appears  to  fall 
within  the  monazite  group.  A  straight  line  fitted  through 
three  uraninite  samples  intersects  the  concordia  curve 
at  1750±30  m.y.  and  368±100  m.y.  and  passes  through  the 
most  reversely  discordant  monazite  (23M).  The  age 
obtained  from  the  upper  intersection  is  the  more 
reliable  one  and  shows  the  "true"  age  of  the  suite 
of  samples.  In  this  case  the  upper  intersection  age  is 
controlled  by  the  concordant  uraninite  point  and  is  in 
good  agreement  with  the  age  obtained  from  the  whole  rock 
lead-lead  data  (to  be  discussed  later) .  This  could  mean 
that  1750  m.y.  represents  the  completion  time  of  the 
uranium  mineralization  which  is  marked  by  the  closing  of 
the  uraninite  crystalline  system.  An  age  of  368  m.y. 
obtained  from  the  lower  intersection  with  the  concordia 
curve,  suggesting  an  episodic  event  that  caused  the  lead 
loss,  is  not  supported  by  any  geological  or  geochrono- 
logical  evidence  up  to  this  time.  Therefore  it  is  assumed 
that  this  age  is  artificial  and  a  probable  result  of  the 
continuous  lead  diffusion  from  the  mineral  lattice.  In 
fact,  the  uraninite  data  define  a  straight  line  rather 
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poorly  (MSWD=19 . 3) .  Comparing  this  data  to  the  contin¬ 
uous  diffusion  curve  with  the  the  constant  coefficient  of 
diffusion:  D(t)=Do  (Wasserberg  1963),  for  a  primary  age 

of  1750  m.y.,  it  becomes  apparent  that  three  uraninite 
points  define  the  continuous  diffusion  curve  much  better 
than  the  straight  line.  Unfortunately,  since  no  more 
uraninite  data  are  available,  a  statement  suggesting  that 
these  uraninites  might  represent  a  good  case  of  continuous 
diffusion  lead  loss  would  be  rather  presumptuous. 

Monazite-r ich  fractions  on  the  other  hand  form  an 

almost  perfect  line  (MSWD=1.1).  Leached  uraninite  (20U) 

has  been  excluded  from  the  straight  line  calculation, 

although  it  fits  very  well  with  this  data  and  does  not 

alter  the  results.  The  line  intersects  the  concordia 

curve  at  1797±3  m.y.  and  at  55±40  m.y.  The  immediate 

implication  of  this  distribution  is  that  U-Pb  systems  in 

monazite-rich  fractions  became  closed  at  1797  m.y.  which 

is  distinctly  different  from  the  cogenetic  uraninites. 

The  lower  intersection  at  55±40  m.y.  suggests  a  recent 

episodic  lead-uranium  loss  or  one  of  the  theoretically 

limiting  cases  (Wetherill  1963)  of  continuous  diffusions 

where  the  intrinsic  diffusion  coefficient  D  is  very 

o 

small  and  the  radiation  damage  of ■ the  lattice  is  somehow 
"frozen  in",  so  that  defects  will  have  negligible  mobility. 
In  that  way,  the  diffusion  coefficient  will  not  increase 
until  the  system  is  raised  to  a  somewhat  higher  tempera¬ 
ture  or  lower  pressure,  sometime  in  the  recent  past,  when 
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the  mobility  of  the  defects  becomes  significant  and 
therefore  the  diffusion  coefficient  is  increased.  The 
result  of  this  quasi-episodic  loss  would  be  the  continuous 
diffusion  curve  that  approaches  the  line  of  recent  episodi 
lead/uranium  loss.  Whatever  the  reason  for  the  lower 
intersection,  the  upper  intersection  should  reflect  the 
"true"  age  of  these  minerals.  Referring  back  to  the 
statement  of  Lewrey  et  al.  (1978)  regarding  the  duration 
of  the  high  grade  metamorphic  condition  in  the  Cree  Lake 
zone  (see  page  7  of  this  thesis),  two  different  ages 
obtained  for  oogenetic  uraninites  and  monazite-rich 
phases  suggest  that  parts  of  the  lead-uranium  system  of 
the  Charlebois  Lake  area  began  closing  at  least  1797  m.y. 
ago  with  crystallization  of  monazites  and  possibly  other 
high  grade  minerals.  The  end  of  the  Hudsonian  orogeny 
in  the  area  was  characterized  by  the  formation  of  urani¬ 
nites,  and  probably  other  late  minerals,  around  1750  m.y. 
ago . 

There  are  other  possible  causes  for  the  observed 
distribution  of  data  points  in  the  concordia  diagram 
in  Figure  5.  These  other  causes  are  suggested  by  the 
experimental  studies  on  the  mobility  of  lead,  uranium 
and  thorium  in  monazites  and  other  uranium  minerals 
(Komlev  et  al .  1964,  Starik  and  Lazarev  1964,  Burger  et 
al.  1967,  Shestakov  1972,  Cumming  and  Rimsaite  1979),  and 
the  theoretical  basis  for  diffusion  of  lead  and  uranium 
provided  by  Wetherill  (1963)  and  Wasserburg  (1963) . 
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Wetherill  (1963)  and  Wasserburg  (1963)  have  proposed 
that  distinctly  different  patterns  for  pure  continuous 
diffusional  losses  of  lead  or  uranium  are  to  be  expected 
in  the  concordia  diagram,  where  the  slope  of  the  pure  lead 
diffusion  line  is  always  greater  than  that  for  pure  uranium 
diffusion.  Mixed  lead  and  uranium  diffusion,  where  the 
diffusion  coefficients  for  lead  and  uranium  are  propor¬ 
tional,  will  produce  a  line  of  intermediate  slope.  In  a 
limiting  case,  if  the  diffusion  coefficient  for  the  parent 
isotope  is  equivalent  to  the  diffusion  coefficient  for  the 
daughter  isotope,  the  resulting  data  point  would  be  con¬ 
cordant,  i.e.  plotting  on  the  concordia  curve. 

Burger  et  al.  (1967)  and  Korolev  et  al.  (1964)  have 
demonstrated  that  various  amounts  of  lead/  uranium  and 
even  thorium  can  be  leached  from  monazite  under  very  mild 
laboratory  conditions.  This  seems  to  be  supported,  to 
some  extent,  in  this  study  by  the  result  of  the  uraninite 
sample  (20U)  which  was  separated  from  the  rest  of  the 
heavy  minerals  by  leaching.  This  sample  appears  to  have 
gained  some  lead  (possibly  a  minor  amount  of  uranium,  too) 
during  the  leaching  process.  Its  position  in  the  concordia 
diagram  depends  on  the  relative  amounts  of  lead  and  uranium 
leached  from  the  monazite  and  other  admixed  minerals,  and 
on  the  total  amount  of  uraninite  present.  It  appears  that 
the  amount  of  leached  lead  added  to  the  uraninite  was 
larger  than  the  amount  of  leached  uranium,  since  the 
addition  of  lead  forces  the  uraninite  towards  the  uranium 
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loss/lead  gain  field,  making  it  reversely  discordant.  At 
the  same  time  the  monazite-rich  fraction  (20M)  of  this 
sample  would  move  from  its  original,  presumably  reversely 
discordant  position  along  the  modern  episodic  lead  loss 
line.  The  distance  of  movement  would  be  determined  by 
the  relative  proportions  of  leached  lead  and  uranium. 

The  other  monazite-rich  fractions  were  washed  in  warm 
2N  HNO^  for  one  hour  before  analysis,  without  immediate 
realization  that  some  lead,  uranium  and  possibly  thorium 
might  be  lost  in  the  process.  If  an  assumption  is  made 
that  various  amounts  of  lead,  uranium  and  thorium  have  been 
indeed  leached  from  the  monazite-rich  fractions  during  the 
washing  process,  then  their  presently  observed  distribution 
is  artificial.  Assuming  they  were  originally  very  closely 
spaced  in  the  vicinity  of  the  most  discordant  monazite 
fraction  (23M),  being  along  the  lead-uranium  continuous 
diffusion  line  represented  by  three  uraninite  samples,  the 
loss  of  lead  and  uranium  would  alter  their  position  and 
move  them  along  the  modern  episodic  lead-uranium  loss  line, 
possibly  achieving  the  observed  distribution. 

This  suggestion  appears  to  be  rather  unlikely  because 
it  assumes  that  all  monazite-rich  fractions  had  almost 
the  same  lead-uranium  composition  as  23M.  The  only  plaus¬ 
ible  fact  about  it  is  that  it  does  not  require  monazite- 
rich  fractions  to  behave  as  closed  systems  under  still 
relatively  high  pressure/temperature  conditions  for  almost 
fifty  million  years  while  U-Pb  systems  of  their  cogenetic 
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uraninites  were  open. 

The  thorium-lead  ages  are  discordant  and  do  not 
substantiate  either  possibility.  It  is  clear  that  addi¬ 
tional  work  of  this  kind  is  needed  to  definitely  establish 
the  age  relationship  of  these  cogenetic  minerals. 

Lead-Lead  Results 

The  analytical  results  obtained  on  the  whole  rock 
samples  and  six  feldspar  separates  are  listed  in  Table  V 
and  the  data  are  plotted  and  presented  in  Figures  6-13. 

All  of  the  data  have  been  corrected  for  fractionation. 

Two  sigma  errors  of  0.7°/oo,  l°/ot>  and  1.4°/00  have  been 
assigned  to  2 0 6 Pb/ 2 0 4 Pb ,  207Pb/204Pb  and  208Pb/204Pb 

ratios,  respectively,  based  on  the  repeated  measurements 
of  NBS  SRM-981  common  lead  standard.  A  detailed  discussion 
of  error  estimates,  analytical  procedures  and  fractiona¬ 
tion  corrections  is  given  in  Appendix  II.  On  all  of  the 
diagrams,  Stacey  and  Kramers  (1975)  growth  curve  is  drawn 
as  a  reference  growth  curve.  In  some  of  the  diagrams  the 

elipses  are  used  for  presenting  data  points.  Their  size 
represents  the  error  limits  at  three  sigma  level. 

A  distribution  of  the  206Pb/204Pb,  207Pb/204Pb  and 

208Pb/204Pb  ratios  obtained  from  the  samples  of  the  Pegasus 
Lake  area  is  displayed  in  Figures  6  and  11.  In  both  of 
these  standard  diagrams  the  data  exhibits  wide  scatter 
which  is  far  beyond  the  assigned  errors.  Some  relatively 
weak  correlation  between  the  isotopic  ratios  and  geo¬ 
graphical  and  geological  locations  appears  to  be  evident 
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in  Figure  6.  Analytical  data  points  obtained  on  the 
samples  from  the  central  areas  of  the  Pegasus  Lake  anti¬ 
cline  seem  to  generally  plot  to  the  left  on  the  diagram. 

The  radiogenic  character  of  the  data  points  increases 
towards  the  shore  of  the  Dramnitzke  Bay,  i.e.  outer 
margins  of  the  anticline.  Samples  of  the  mineralized 
granodioritic  granofels  (18,19,20),  with  206Pb/204Pb 
ratios  of  up  to  240  plot  outside  this  diagram  and  lie  on 
line  V  whose  extension  is  shown  as  a  broken  line  in  this 
diagram . 

The  somewhat  fan-like  distribution  of  data  in  Figure 
6  is  indicative  of  a  complex  history  for  these  rocks.  Lead- 
uranium  systems  from  the  rocks  with  a  simple  history  (i.e. 
rocks  which  were  emplaced  sometime  in  the  past  and  have 
not  been  disturbed  until  the  present)  exhibit  colinear 
behaviour  in  the  standard  lead-lead  diagrams.  Systems 
with  a  more  complex  history  occasionally  exhibit  straight 
line  patterns  and  the  slopes  of  these  straight  lines  might 
often  reflect  the  time  of  the  last  metamorphic  event  but 
that  is  only  true  in  a  few  special  cases  (discussed  in 
Appendix  I) .  Colinear  behaviour  is  observed  only  when  the 
changes  are  brief,  episodic  and  complete,  i.e.  when  the 
system  exhibits  closed  system  behaviour  and  when  all  the 
chemical  and  isotopic  equilibration  reactions  are  completed 
during  a  brief  metamorphic  event.  If  the  system  remains 
open  for  partial  uranium  or  lead  migration  across  the 
system's  boundaries,  colinearity  is  seldom  observed. 
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Figure  6.  Pb-Pb  diagram  for  the  whole  rock  samples  from  the  "Pegasus  Lake  Area"  For 

identification  of  the  specimens  refer  to  Figure  4  and  Table  V.  Data  ellipses 
are  drawn  at  3a  error  level. 


45 


Providing  a  sufficient  number  of  data  is  available, 
some  relatively  useful  information  regarding  the  history 
of  the  system  may  be  derived  from  the  scatter  diagram. 

The  slope  age  calculated  from  the  line  representing  the 
lower  boundary  of  the  scatter  often  approaches  the  time 
of  the  last  metamorphic  event.  The  distribution  of  data 
points  around  the  apex  of  the  "fan"  is  occasionally 
indicative  of  the  time  of  the  original  formation. 

The  lower  boundary  of  the  scatter  in  the  207Pb/204Pb 
vs.  206Pb/204Pb  diagram  (Figure  6)  is  represented  by  a 
line  fitted  through  seven  points  with  the  slope  of  0.10855 
and  MSWD  of  2.45.  The  corresponding  slope  age  is  1775180 
m.y.  at  two  sigma  error  level.  For  reference  this  line  is 
termed  "Line  I".  The  intersection  of  this  line  with 
Stacey  and  Kramers  growth  curve  at  about  2100  m.y.  and 
distribution  of  the  data  points  suggests  that  the  original 
emplacement  took  place  sometime  before  2100  m.y.  ago  and 
subsequent  remobilization  occurred  at  around  1775  m.y.  ago. 

It  will  only  be  mentioned  now  (and  discussed  later) 
that  three  samples  of  granodioritic  granofels  (18,19,20) 
and  two  samples  from  the  "Old  Camp  Area"  plot  on  line  V 
which  represents  the  upper  limit  of  the  observed  distribu¬ 
tion  in  Figures  6  and  7.  The  slope  of  this  line  yields 
an  age  of  1772±8  m.y.  (Figure  8) . 

Isotopic  results  obtained  on  the  gneissic  samples 
from  the  "Old  Camp  Area"  are  plotted  in  Figures  7  and  11. 
The  data  in  the  207Pb/204Pb  vs.  206Pb/204Pb  diagram  is 
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distributed  along  three  parallel  trend  lines  (with  the 
exception  of  sample  27).  Straight  line  fittings  of  these 
data  yield: 


Line  II:  Slope=0 . 10 6 7 8 
Line  III:  Slope=0 . 10 72 7 
Line  IV:  Slope=0 . 10708 


MSWD=  7.91  T  =1745179  m.y.  (2o) 
s 

MSWD=2 . 55  T  =1754126  m.y.  (2a) 
s 

Two  point  line  Tg=1751  m.y. 


As  in  the  case  of  the  "Pegasus  Lake  Area"  some 
correlation  between  the  distribution  of  the  data  in  Figure 
7  and  the  geographic-geological  location  of  the  samples 
can  be  made.  Most  radiogenic  samples  were  collected  along 
the  shore  of  the  Dramnitzke  Bay.  The  radiogenic  character 
of  the  specimens  decreases  towards  the  apparent  fault  and 
amphibolite  area  (samples  26  and  27).  Northwest  of  the 
apparent  fault  the  correlation  is  not  quite  clear,  but  it 
seems  that  the  rocks  adjacent  to  the  fault  area  are  more 
radiogenic  and  the  radiogenic  character  decreases  towards 
the  more  felsic  granitic  gneiss  (samples  32,33,34).  Sampl 
34  probably  owes  its  position  on  the  207Pb/204Pb  vs. 
206Pb/204Pb  diagram  to  the  radiogenic  development  of  its 
U-Pb  system  before  the  last  metamorphic  event.  In  order 
to  have  such  a  high  207Pb/204Pb  ratio  and  relatively  low 
206Pb/204Pb  ratio,  the  y  value  of  its  Pb-U  system  before 
the  last  metamorphic  event  must  have  been  higher  than  9.74 
(which  is  Stacey  and  Kramers  value  for  the  earth's  y2 
environment  (Stacey  and  Kramers  1975))  and  it  must  have 
lost  most  of  its  uranium  during  the  last  metamorphic  event 


Therefore,  the  y  value  for  this  system  was  low  in  the 
last  stage  of  its  history. 

Samples  22  and  34  form  line  IV  which  represents  the 
upper  limit  to  the  distribution  of  the  data  from  the  "Old 
Camp  Area".  When  these  data  are  combined  with  the  grano- 
dioritic  granofels  data  (samples  18,19,20)  and  plotted 
together  (Figure  8)  the  linear  trend  with  the  slope  of 
0.10840  is  obtained  (line  V).  Although  the  MSWD  of  12.58 
suggests  that  the  data  do  not  determine  the  line  very  well, 
the  calculated  slope  age  at  two  sigma,  taking  into  account 
the  error  due  to  the  imperfect  fit,  yields  an  age  of 
1772  ±8  m . y . 

Three  parallel  line  distributions  of  data  points,  the 
lowest  line  intersecting  the  reference  growth  curve  at 
around  2200  m.y.,  strongly  suggest  that  these  rocks  might 
have  developed  according  to  the  three  stage  evolution  model 
with  a  constant  f^  fractionation  factor  (Gale  and  Mussett , 
1973) .  The  position  of  the  lowest  line  indicates  that  the 
original  emplacement  of  these  rocks  must  have  happened 
more  than  2200  m.y.  ago  and  that  subsequent  remobilization 
took  place  between  1700  and  1800  m.y.  ago.  It  appears  that 
the  Pb-U  system  of  the  "Old  Camp  Area"  was  fractionated 
among  at  least  three  major  subsystems  at  the  time  of  the 
original  emplacement  and  then  it  was  refractionated  during 
the  Hudsonian  orogeny,  producing  the  presently  observed 
distribution  of  three  parallel  lines. 

Keeping  in  mind  the  conclusions  drawn  by  Morra  (1977) 
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Figure  8.  Pb-Pb  diagram  for  the  whole  rock  granodiori tic  granofels  (circles)  and  some 
migmatite  (squares)  members  of  the  Charlebois  Lake  Complex.  For  identifica¬ 
tion  of  the  specimens  refer  to  Figure  4  and  Table  V. 


regarding  the  integrity  of  the  Charlebois  Lake  complex, 
some  similarities  to  the  results  from  the  "Pegasus  Lake 
Area"  and  the  "Old  Camp  Area"  are  obvious.  The  slopes  and 
slope  ages  obtained  from  line  1  and  line  2  (1775±80  m.y. 
and  1745±79  m.y.)  and  the  intercepts  of  these  lines  with 
reference  growth  curve  (2100  m.y.  and  2192  m.y.)  are  well 
within  the  assigned  error  and  in  both  cases  the  lines 
represent  the  lower  limits  of  the  observed  distribution 
in  the  207Pb/204Pb  vs.  206Pb/204Pb  diagrams.  The  upper 
limit  of  the  distribution  in  both  cases  is  represented  by 
line  V  with  the  slope  age  of  1772±8  m.y.  This  effectively 
puts  the  upper  and  lower  limits  on  the  \i  environment  of 
the  Charlebois  Lake  complex  for  the  time  period  between 
the  present  and  the  Hudsonian  metamorphic  event.  The 
geologic  evidence  and  above  distribution  make  it  reason¬ 
able  to  assume  that  Pb-U  systems,  from  the  "Pegasus  Lake 
Area"  and  the  "Old  Camp  Area"  had  also  a  common  history 
before  the  Hudsonian  time. 

In  an  attempt  to  gain  some  information  about  the 
time  of  the  original  emplacement  of  the  Charlebois  Lake 
complex  and  the  history  of  the  intermediate  stage, 
feldspars  from  six  of  the  least  radiogenic  samples  were 
separated  and  analysed.  The  chosen  samples  were:  three 
specimens  from  the  "Pegasus  Lake  Area"  (1,  8  and  13)  and 
three  from  the  "Old  Camp  Area"  (27,  29  and  34) . 

Under  ideal  conditions  lead  isotopic  composition  in 
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the  feldspars  should  be  representative  of  the  U  environ¬ 
ment  in  the  respective  Pb-U  systems  at  the  beginning  of 
the  last  evolutionary  stage.  Low  affinity  of  the  feldspar 
lattice  for  uranium  means  that  the  isotopic  composition 
of  lead  in  the  system  at  time  i.e.  the  Hudsonian 

metamorphic  event,  would  be  "frozen  in"  and  no  in  situ 
generation  of  radiogenic  lead  would  be  possible.  That 
means  that  feldspars  behave  to  some  extent  as  galena 
(Sinha  1969),  i.e.  the  lead  evolution  should  stop  at  time 
and  the  data  points  should  form  a  linear  trend  in  the 
2Q7Pb/204Pb  vs.  206Pb/204Pb  diagram.  The  slope  of  this 
line  would  be  determined  by  the  time  of  the  original 
emplacement  of  the  rocks  (t^)  and  subsequent  metamorphic 
event  ( ) • 

The  slopes  of  the  lines  connecting  the  data  points 
of  the  whole  rocks  to  that  corresponding  feldspars  should 
yield  the  age  of  the  metamorphic  event. 

Lead  isotopic  ratios  of  feldspar-whole  rock  pairs 
from  the  Charlebois  Lake  complex  are  plotted  in  the 
207Pb/204Pb  vs.  206Pb/204Pb  diagram  (Figure  9).  All 
feldspar  separates  were  washed  in  warm  4N  HNO^  for  two 
hours  before  the  lead  was  separated  (described  in 
Appendix  II).  Samples  treated  in  this  manner  have  the 
subscript  "F"  added  to  the  sample  I.D.  number.  Later 
on,  during  this  study,  the  remaining  amount  of  one  of  the 
feldspar  separates  (13F)  was  leached  with  a  warm  mixture 
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consisting  of  20  ml  concentrated  HNO^  and  0.5  ml  con¬ 
centrated  HF  for  one  hour.  The  residue  (13FR)  was  then 
analysed  for  lead  isotopes.  This  data,  along  with  the 
data  from  the  other  feldspars,  is  presented  in  Table  V. 

As  expected,  all  feldspar  ratios  are  less  radiogenic 
than  their  whole  rock  counterparts  but  they  exhibit  large 
scatter  instead  of  forming  a  linear  trend.  The  lines 
connecting  the  pairs  are  far  from  being  parallel  and  the 
ages  calculated  from  these  lines  range  from  around 
1500  m.y.  to  2500  m.y.  Even  though  the  errors  associated 
with  the  slopes  and  therefore  with  the  calculated  ages 
are  largely  due  to  the  small  spread  in  the  feldspar-whole 
rock  206Pb/204Pb  ratios,  this  cannot  be  considered  as  the 
only  reason  for  such  a  large  variation  in  the  calculated 
slope  ages. 

Earlier  isotopic  studies  of  old  Precambrian  feldspars 
have  shown  that  feldspars  could  be  useful  in  whole  rock 
lead  dating  (Davis  1978,  Vidal  et  al.  1980)  but  they  often 
yield  younger  model  ages  (Sinha  1969,  Catinzaro  and  Gast 
1960,  Doe  et  al .  1965) .  The  lead  isotopic  composition 
of  feldspars  commonly  approaches  that  of  associated  lead 
ores,  but  it  is  generally  more  radiogenic.  Doe  et  al. 
(1965)  suggested  that  gain  or  exchange  of  feldspar  lead 
with  other  rock  or  mineral  lead  or  admixing  of  lead  from 
fluids  are  possible  mechanisms  causing  the  increased 
radiogenic  character  of  the  feldspar  leads. 
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Sinha  (1969)  in  the  lead  volatilization  study  of 
feldspars,  indicated  that  radiogenic  lead  can  be  easily 
removed  from  feldspar  by  preferential  volatilization 
since  this  lead  is  concentrated  along  the  grain  boundaries 
or  in  distorted  lattice  sites.  Ludwig  and  Silver  (1977) 
recognized  two  types  of  unsupported  radiogenic  lead 
present  in  feldspars.  One  has  normal  isotopic  composi¬ 
tion  with  relative  abundances  of  206Pb,  207Pb  and  208Pb 
consistent  with  the  age  and  uranium- thorium  content  of  the 
rock.  The  other  type  of  unsupported  radiogenic  lead  is 
apparently  pure  206Pb  isotope,  derived  from  long-term 
migration  and  accumulation  of  the  radioactive  daughters 
of  238U,  namely  222Rn  and/or  226Ra. 

If  it  is  accepted  that  various  amounts  of  unsupported 
lead  are  present  in  the  analysed  feldspars  from  the  Charle- 
bois  Lake  complex,  then  it  is  clear  that  lead  data  will 
not  form  a  linear  trend  but  scattered  points.  The  presence 
of  pure  206Pb  isotope  component  would  certainly  affect  the 
position  of  the  feldspar  data  points  in  the  207Pb/204Pb  vs. 
206Pb/204Pb  diagram  and  the  slopes  of  the  lines  connecting 
these  points  and  their  corresponding  whole  rock  data 
points,  the  result  being  the  horizontal  shift  in  the 
feldspar  data  and  higher  slopes. 

Isotopic  ratios  obtained  from  the  residue  (13FR)  of 
the  leached  feldspar  (13F)  seem  to  support  the  above 
argument,  since  the  isotopic  results  of  the  residue  (13FR) 
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are  less  radiogenic  than  those  of  the  total  feldspar 
separate  (13F)  and  the  shift  in  the  206Pb/204Pb  appears 
to  be  horizontal.  Unfortunately  there  is  not  enough 
data  of  this  kind  available  at  present  to  support  this 
finding.  Therefore,  it  is  immaterial  whether  the 
observed  feldspar  distribution  in  Figure  9  is  a  result 
of  unsupported  radiogenic  lead  or  incomplete  equilibra¬ 
tion  and  homogenization  of  the  whole  rock  Pb-U  system 
with  the  feldspar  Pb-U  system  during  the  metamorphic 
event.  Firm  conclusions  regarding  the  history  of  the 
Charlebois  Lake  complex  before  the  Hudsonian  metamorphic 
event  cannot  be  drawn  from  the  data  presented  so  far. 

Nevertheless,  with  the  few  reasonable  assumptions 
regarding  the  feldspars,  a  very  plausible  model  history 
of  the  Pb-U  systems  in  the  Charlebois  Lake  area  can  be 
obtained.  It  is  quite  reasonable  to  assume  that  the 
feldspars  least  affected  by  the  possible  presence  of 
unsupported  radiogenic  lead,  i.e.  the  ones  closest  to 
the  ideal  behaviour,  would  plot  to  the  far  left  on  the 
207Pb/204Pb  vs.  206Pb/204Pb  diagram.  A  slope  of  0.27401 
(MSWD-8.95)  is  obtained  when  the  line  is  fitted  through 
four  feldspar  points  (Figure  9).  Since  most  of  the 
unsupported  radiogenic  lead  is  presumably  removed  from 
the  feldspar  13FR  and  the  other  three  feldspars  might 
contain  some  of  this  lead,  the  calculated  slope  very 
likely  represents  the  maximum  possible  slope.  The  slope 
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of  the  line  is  determined  by  the  time  of  the  original 
emplacement  (t^)  and  the  time  of  the  subsequent  meta- 
morphic  event  (t^)-  Therefore  for  the  assigned 
metamorphic  age  =  1775  m.y.,  the  maximum  initial  age 
t^  of  2470  m.y.  is  obtained.  Knowing  the  times  t^  and 
t^  and  the  slope  of  the  feldspar  line,  the  average  y^ 
value  is  calculated  to  be  9.85. 

A  graphical  representation  of  the  postulated  three 
stage  model  history  with  f^  =  const,  of  the  Charlebois 
Lake  complex  is  given  in  Figure  10. 

The  overall  Pb-Pb  data  suggests  that  during  the 
period  before  the  emplacement  of  the  Charlebois  Lake 
complex,  i.e.  before  2470  m.y.,  uranium  and  lead  of  this 
complex  were  developing  in  an  environment  with  average 
y^  =  9.85.  As  early  as  2470  m.y.  ago  y^  environment 
was  rehomogenized  and  refractionated  by  erosion, deposi¬ 
tion,  and  sedimentation  into  at  least  three  second-stage 
y2  environments.  A  Hudsonian  metamorphic  event,  nominally 
taken  here  to  have  occurred  1775  m.y.  ago,  initiated  the 
third  stage  of  Pb-U  evolution  and  was  responsible  for  the 
presently  observed  distribution  of  data  points  in  the 
207Pb/204Pb  vs.  206Pb/204Pb  diagram  (Figure  10).  Devia¬ 
tion  from  the  ideal  three-stage,  f^  =  const,  model 
history,  expressed  by  the  imperfect  line  fit  of  the  data 
and  general  scatter  could  be  attributed  to  the  open  system 
behaviour.  It  was  suggested  earlier,  by  the  uraninite 
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and  monazite  U-Pb  data,  that  the  global  Pb~U  system  in 
the  Charlebois  Lake  area  had  remained  partially  open  for 
approximately  fifty  million  years  during  the  Hudsonian 
orogeny.  This  appears  to  be  particularly  evident  in  the 
Pegasus  Lake  area  where  the  data  points  generally  scatter 
between  the  limiting  line  I  and  line  V.  It  is  evident 
that  the  most  radiogenic  samples,  i.e.  granodior itic 
granofels  and  migmatites  which  plot  on  line  V  have 
partially  evolved  from  the  high  environment  of  the 
intermediate  stage.  That  might  suggest  that  the  uranium 
concentration  in  the  original  sediments  was  highest  near 
the  boundary  between  arkosic  and  more  mafic  calcareous 
and  tufaceous  sediments,  and  lowest  at  the  bottom  of  the 
sedimentary  sequence  which  is  seen  now  as  the  core  of  the 
anticlines  in  the  area.  In  addition,  not  very  large 
scale  but  rather  local  migration  and  concentration  of 
uranium  during  the  Hudsonian  orogeny  has  produced  uranium 
mineralization  which  is  associated  with  granodior itic 
granofels  and  migmatites.  This  seems  to  be  supported  by 
the  linearity  of  the  data  trends  in  Figure  10.  If  large- 
scale  uranium  remobilization  had  taken  place  during  the 
Hudsonian  time,  the  colinear  behaviour  of  the  data  would 
have  been  destroyed.  Partial  mixing  of  the  lead  and 
uranium  from  different  subsystems  during  the  migration 
would  have  produced  general  scatter.  That  might  well  be 
the  case  with  some  of  the  granitic  gneiss  samples  from 
the  Pegasus  Lake  anticline. 
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Samples  from  the  center  of  the  anticline  which 
evolved  in  the  relatively  low  y 2  environment  in  the 
intermediate  stage,  seem  to  have  lost  most  of  their 
uranium  and  have  also  undergone  most  complete  re¬ 
equilibration  during  the  last  metamorphic  event.  As  a 
result  they  form  the  lower  boundary  of  the  distribution. 
Samples  from  the  outer  margins  of  the  anticline  have 
probably  been  involved  in  the  partial  migration  of  lead 
and  uranium  which  has  resulted  in  general  scatter  of 
the  data  from  that  area  since  rehomogenization  has  not 
been  completed. 

Figure  11  shows  the  208Pb/204Pb  versus  206Pb/204Pb 
ratios  of  the  Charlebois  Lake  complex.  The  data  show 
wide  scatter  as  is  often  the  case  in  trace  lead  studies. 
In  general  the  scatter  of  the  data  points  indicates 
considerable  variation  in  the  Th-U  ratio  and  therefore 
it  is  consistent  with  the  proposition  made  above  that 
several  events  have  affected  the  development  of  the  lead 
in  the  rocks  of  the  Charlebois  Lake  complex. 

Late  Felsic  Dyke 

A  fifteen  foot  wide  dyke  which  crosscuts  the  rocks 
in  the  "Old  Camp  Area"  and  grades  into  the  small,  very 
coarse  grained  pink  pegmatitic  irregular  body  on  the 
north  side  of  the  bay  (Figure  4) ,  was  sampled  at  six 
locations  (35  to  40).  Samples  35  and  36  were  collected 
one  and  three  feet,  respectively,  from  the  contact  with 
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206  7  204 
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Figure  11.  Pb-Pb  diagram  for  the  whole  rock  samples  of 
various  rock  units  from  the  Charlebois  Lake 
Complex . 
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the  country  rock.  Sample  37  is  from  the  center  of  the 
dyke  and  the  other  three  samples  were  collected  elsewhere 
from  the  dyke. 

Lead  analytical  data  are  presented  in  Table  V  and  in 
Figures  12  and  13.  Data  plotted  in  the  207Pb/204Pb  versus 
206Pb/204Pb  diagram  display  a  linear  trend  with  the  range 
of  206Pb/204Pb  ratios  between  17  and  19.5.  When  a 
straight  line  is  fitted  to  the  data  it  yields  a  slope  of 
0.08152  and  an  MSWD  of  3.39.  The  corresponding  slope  age 
is  1235±320  m.y.  The  line  intersects  Stacey  and  Kramers 
reference  growth  curve  at  -8±230  m.y.  and  1239±300  m.y. 
at  a  two  sigma  error  level.  Due  to  the  small  spread  in 
the  206Pb/204Pb  ratios,  the  small  number  of  samples  and 
general  scatter,  the  data  points  define  the  line  rather 
poorly.  As  a  result  the  calculated  slope,  slope  age  and 
intercepts  with  the  growth  curve  are  associated  with 
fairly  large  error  and  the  interpretation  of  this  data 
becomes  only  tentative. 

It  is  interesting  to  note  that  the  line  intersec¬ 
tions  with  Stacey  and  Kramers  reference  growth  curve 
indicate  that  the  source  of  this  dyke  has  a  y  value  very 
similar  to  the  y  value  of  the  upper  mantle  as  suggested 
by  Stacey  and  Kramers  (1975) .  Unfortunately  this 
suggestion  is  not  supported  by  the  distribution  of  the 
data  points  in  the  208Pb/204Pb  versus  206Pb/204Pb 
diagram  (Figure  13)  nor  by  Rb-Sr  data  which  will  be 
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Figure  12.  Pb-Pb  diagram  for  the  whole  rock  samples  of  the  pink  felsic  dyke.  For 
identification  of  the  specimens  refer  to  Figure  4  and  Table  V.  Data 
ellipses  are  drawn  at  3a  error  level. 
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Figure  13. 

206  ,  204 

Pb/  .  Pb 

Pb-Pb  diagram  for  the  whole  rock  samples 
of  the  pink  felsic  dyke. 
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discussed  later  in  this  chapter.  What  might  be  said  is 
that  the  source  of  the  felsic  dykes  and  other  irregular 
pegmatitic  bodies  is  probably  remelted  rocks  of  the 
lower  crust.  The  isotopic  data  confirm  the  geologic 
observation  that  emplacement  of  this  and  other  similar 
dykes  and  pink  pegmatitic  irregular  bodies  has  happened 
in  the  post  metamorphic  time,  possibly  as  late  as  around 
1200  m.y.  ago. 

The  period  of  time  between  1500  to  around  1100  m.y.  ago 
is  tectonically  a  relatively  quiescent  and  economically 
important  period  in  the  history  of  the  Churchill  province. 

It  should  be  pointed  out  that  during  this  time  major 
uranium  deposits  of  northern  Saskatchewan  such  as  Rabbit 
Lake,  Midwest  Lake,  Key  Lake  and  others  have  been 
emplaced  (Cumming  and  Rimsaite  1979,  Cumming,  Baadsgaard, 
Hoeve  and  Morton  (in  preparation),  and  Cumming,  Worden 
and  Baadsgaard  (in  preparation)).  Some  of  the  diabase 
dykes  of  Churchill  province  have  also  been  emplaced  during 
this  period  of  time  (Fahrig  and  Wanless  1963,  Burwash  et 
al.  1962).  A  detailed,  large  scale  geochronological  study 
of  these  and  felsic  dykes  from  other  areas  in  the  Churchill 
province,  along  with  the  diabase  dykes,  would  probably 
prove  to  be  very  beneficial  for  a  clearer  understanding 
of  the  Churchill  province  history  between  1500  m.y.  and 


1100  m.y.  ago. 
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Rubidium-Strontium  Results 

The  analytical  results  obtained  on  the  whole  rock 
samples  are  listed  in  Table  VI  and  the  data  are  plotted 
in  Figures  14-16. 

Two  sigma  errors  of  0.3°/oo  and  1%  have  been  assigned 
to  87Sr/86Sr  and  87Rb/86Sr  ratios,  respectively,  based  on 
the  repeated  measurements  of  NBS  SRM-987  common  strontium 
standard  and  the  generally  accepted  accuracy  level  in  the 
laboratory.  A  detailed  discussion  of  error  estimates  and 
analytical  procedures  is  given  in  Appendix  II.  In  some 
of  the  diagrams  the  elipses  are  used  for  presenting  the 
data  points,  their  size  represents  the  data  with  their 
respective  errors  at  four  sigma  level. 

For  ease  of  interpretation  the  data  are  again 
considered  in  three  groups  as  was  the  case  with  the 
lead-lead  data,  i.e.  "Pegasus  Lake  Area",  "Old  Camp  Area" 
and  the  late  felsic  dyke. 

All  of  the  rock  samples  analysed  for  lead  isotopes 
were  also  used  for  Rb-Sr  analysis,  in  an  attempt  to  shed 
some  more  light  on  the  history  and  the  nature  of  the 
Charlebois  Lake  complex  and  its  precursors.  It  was  hoped 
that  the  whole  rock  Rb-Sr  analysis  would  provide  addi¬ 
tional  evidence  for  the  age  of  the  precursors  and  there¬ 
fore  substantiate  the  claims  made  on  the  basis  of  the 
lead-lead  data.  Unfortunately  the  data  presented  here 
exhibits  even  greater  scatter  than  the  lead-lead  data. 
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The  Hudsonian  metamorphic  event  appears  to  have  been 
strong  enough  to  wipe  out  the  record  of  the  time  when  these 
rocks  v/ere  originally  emplaced,  but  not  strong  enough  to 
completely  reset  the  whole  rock  Rb-Sr  clock  in  this  area. 

The  Rb-Sr  data  in  all  of  the  standard  isochron 
diagrams  exhibits  wide  scatter  so  that  no  firm  geochrono- 
logical  conclusions  were  possible.  Close  inspection  of 
the  data  in  the  scatter  diagrams  coupled  with  the  evidence 
obtained  from  the  lead  analyses  and  other  geological 
studies  provide  some  indications  of  the  probable  history 
of  the  Rb-Sr  systems  of  these  rocks.  The  purpose  of  the 
discussion  that  follows  is  to  see  what  kind  of  information 
can  be  extracted  from  highly  discordant  Rb-Sr  data  and 
how  closely  it  correlates  with  the  information  obtained 
from  lead-uranium  and  lead-lead  analyses. 

The  data  obtained  on  the  samples  from  the  "Pegasus 
Lake  Area",  including  three  granodior itic  granofels 
specimens  (18,  19  and  20)  are  plotted  in  Figure  14. 

Coarse-grained,  mineralized,  granodior it ic  granofels 
samples  exhibit  a  relatively  wide  spread  in  Rb-Sr  ratios 
lying  along  the  trend  which  would  give  an  initial  87Sr/ 

86Sr  ratio  of  around  0.712  and  slope  age  of  close  to 
1750  m.y.  As  in  the  case  of  lead-lead  data  it  again  forms 
the  upper  limit  to  the  distribution  of  data  in  the  standard 
isochron  plot. 

The  rest  of  the  gneissic  samples  are  fairly  mono- 
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tonous  with  respect  to  the  Rb-Sr  ratios,  exhibiting  a 
limited  range  of  values,  as  would  be  expected  from  their 
mineralogy.  This  is  especially  true  for  samples  from 
the  central  part  of  the  anticline  where  the  gneiss  has  a 
predominantly  gray,  massive  appearance  with  no  foliation 
or  even  lineation.  The  range  of  87Rb/86Sr  values  for 
samples  1  to  13  is  0.23-0.60  (the  points  within  the 
dashed  area  in  Figure  14).  The  rest  of  the  points  which 
do  not  fall  into  this  group  were  collected  from  within 
about  300  meters  from  the  shore  of  Dramnitzke  Bay,  i.e. 
adjacent  to  the  mineralized  granodioritic  granofels  area. 
The  biotite  content  of  these  rocks  is  somewhat  higher 
than  in  the  rocks  from  the  center  of  the  anticline;  the 
colour  is  predominantly  pinkish-gray  and  the  foliation 
is  visible.  The  range  of  87Rb/86Sr  ratios  is  wider  and 
sample  17  has  the  most  radiogenic  value  of  4.2.  A  line 
fitted  through  three  of  these  points  (15,  16  and  17) 
passes  through  the  scatter  displayed  by  the  anticline  core 
samples,  giving  a  slope  age  of  1753  m.y.  and  a  87Sr/86Sr 
initial  ratio  of  0.7066.  It  should  be  noted  that  most  of 
the  other  samples  from  the  scatter  that  seem  to  lie  along 
this  line  are  geographically  close  to  the  granodioritic 
granofels  area. 

An  attempt  to  fit  the  line  through  the  scatter 
displayed  by  the  anticline  core  samples  which  have  a 
small  range  of  87Rb/06Sr  values  and  are  distributed 
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within  the  dashed  area  in  Figure  14  (1-13,  except  #2), 

would  yield  a  slope  age  of  approximately  2400  m.y.  and 
0.701  for  an  initial  87Sr/86Sr  ratio. 

Rb-Sr  systems  for  the  rocks  from  the  "Pegasus 
Lake  Area"  appear  to  show  some  similarities  to  the 
U-Pb  systems.  In  both  cases  the  radiogenic  character 
increases  towards  the  outer  margins  of  the  anticline. 
Granodior i tic  granofels  data  form  the  upper  limit  to 
the  distribution  of  points  in  both  Pb-Pb  and  Rb-Sr 
diagrams.  Rubidium  and  uranium  seem  to  be  enriched 
in  granodiorit ic  granofels,  while  the  rocks  adjacent  to 
the  mineralized  area  exhibit  an  intermediate  character 
between  centrally  located  samples  and  the  samples  from 
the  mineralized  area. 

Data  obtained  from  the  samples  of  the  "Old  Camp 
Area"  clearly  do  not  exhibit  closed  system  behaviour 
either  (Figure  15) .  87Rb/86Sr  values  extend  to  around 

3.5  with  the  majority  of  points  scattering  between  0.0 
and  1.0.  The  data  do  not  display  any  colinearity, 
therefore  any  conclusions  drawn  from  this  data  would 
be  speculative  and  just  an  attempt  to  obtain  confirma¬ 
tion  of  the  rock  history  based  on  the  whole  rock  lead- 
lead  and  uranium-thorium-lead  data. 

Assuming  that  the  mobility  of  the  elements  during 
the  last  metamorphic  event  has  been  limited  to  predomi¬ 
nantly  small-scale  migration  and  exchange,  without  large 
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Figure  15.  Rb-Sr  diagram  for  the  whole  rock  samples  fran  the  "Old  Camp  Area".  "Reference  isochrons 
indicate  the  initial  age  and  metamorphic  age.  For  identification  of  the  specimens  refer 
to  Figure  4  and  Table  VI.  Data  ellipses  are  drawn  at  4a  error  level. 
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contributions  of  rubidium  and  strontium  from  outside, 
a  lower  limit  to  the  87Sr/86Sr  ratio  at  the  end  of  the 
metamorphic  event  (t^)  can  be  established.  This  cannot 
be  lower  than  the  ratio  of  the  system  when  it  was  first 
completely  rehomogenized  (t^).  If  this  value  was  avail¬ 
able,  the  upper  limit  to  the  age  of  the  metamorphic 
event,  which  has  disturbed  the  original  Rb-Sr  system, 
can  be  estimated  by  drawing  a  reference  isochron  that 
starts  at  this  ratio  and  bounds  the  lower  scatter  of 
the  data  points.  In  this  study  the  end  of  the  meta¬ 
morphic  event  is  fairly  well  established  by  the  lead-lead 
and  uranium-lead  results  to  be  around  1750  m.y.  The 
value  sought  here  is  an  87Sr/86Sr  ratio  which  might  have 
been  the  initial  87Sr/86Sr  ratio  during  the  emplacement 
of  the  "Old  Camp  Area"  rocks,  i.e.  the  Charlebois  Lake 
complex.  Fitting  the  line  with  the  slope  age  of  1750  m.y. 
to  the  three  points  representing  the  lower  boundary  of 
the  scatter  in  Figure  20,  an  initial  ratio  of  about  0.701 
was  obtained.  Whatever  the  time  of  the  initial  emplace¬ 
ment  (t^)  of  the  Charlebois  Lake  complex  was,  the  average 
87Sr/86Sr  ratio  could  not  have  been  higher  (or  much  lower 
for  that  matter)  than  0.701.  The  same  value  for  initial 
87Sr/86Sr  ratio  was  obtained  from  the  line  fitted  to  the 
scatter  displayed  by  the  samples  from  the  core  of  the 
Pegasus  Lake  anticline.  Since  this  value,  obtained  in 

8  7  8  6 

two  different  ways,  is  suggestive  of  the  initial  Sr/  Sr 
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ratio  of  precursors  of  the  Charlebois  Lake  complex, 
further  possibilities  for  deducing  the  original  emplace¬ 
ment  age  were  examined. 

For  this  purpose  all  of  the  data  were  considered 
together  and  plotted  in  Figure  16.  Out  of  the  thirty- 
four  samples  analysed,  twenty-seven  had  87Rb/86Sr  ratios 
of  less  than  1.0.  The  upper  and  lower  boundaries  of  the 
data  consist  of  samples  of  granodior itic  granofels  and 
three  samples  from  the  "Old  Camp  Area"  respectively. 
Reference  isochrons  of  1750  m.y.  are  drawn  through  the 
upper  and  lower  boundaries.  This  band-like  distribution 
of  the  Rb-Sr  data  indicates  that  at  the  end  of  the 
Hudsonian  orogeny,  the  global  Rb-Sr  system  of  the 
Charlebois  Lake  complex  was  probably  incompletely  homo¬ 
genized,  with  the  87Sr/86Sr  ratio  ranging  from  around 
0.701  to  0.712  for  the  various  subsystems  of  the  complex. 

It  has  been  already  mentioned  that  metasomatism  in 
the  area  has  not  been  observed  (Morra  1977) .  If  the 
mobility  of  rubidium  and  strontium  within  the  complex  as 
a  whole  was  limited  and  only  local,  then  from  the  mean 
initial  87Sr/86Sr  -  0.705  at  the  end  of  the  last  meta- 
morphic  event  (about  1750  m.y.  ago)  and  the  average 
8  7  Rb / 8  6  S  r  of  0.412  for  the  Charlebois  Lake  complex,  it 
could  be  calculated  that  it  would  take  approximately 
680  m.y.  to  bring  the  87Sr/86Sr  ratio  from  0.701  to  0.705. 
That  effectively  sets  the  emplacement  time  of  the  Charle- 
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R  b  /  Sr 

Figure  16.  Combined  Rb-Sr  diagram  for  the  whole  rock  samples  from  the  "Old  Camp  Area 
and  the  "Pegasus  Lake  Area".  "Reference  isochrons"  indicating  the  age  of 
the  precursors  and  the  metamorphic  age  are  shown. 
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bois  Lake  complex  percursors  at  about  2440  m.y.  ago, 
which  is  in  good  agreement  with  the  age  of  2470  m.y. 
obtained  from  the  feldspar  lead-lead  data  and  the  2400 
m.y.  obtained  from  the  anticline  core  samples. 

The  distribution  of  the  Rb-Sr  data  presented  so  far 
clearly  demonstrates  incomplete  homogenization  during  the 
last  metamorphic  event  or  partial  open  system  behaviour. 
Various  studies  dealing  with  the  remobilization  effects 
of  rubidium  and  strontium  on  mineral  systems,  under 
metamorphic  conditions,  have  been  carried  out  (Baadsgaard 
and  van  Breemen  19  70,  Kesmarky  19  77)  .  It  was  shown  that 
micas  are  readily  susceptible  to  this  kind  of  change. 
Whether  or  not  the  whole  rocks  are  disturbed  will  depend 
upon  the  range  of  migration  of  these  elements.  The 
relative  stability  of  Rb-Sr  systems  in  whole  rocks 
indicates  that  under  many  circumstances  rubidium  and 
strontium  will  become  very  quickly  trapped  in  minerals 
having  low  chemical  potential  for  these  elements  (such 
as  K-feldspar) ,  before  they  have  a  chance  to  migrate  very 
far.  However,  other  studies  have  demonstrated  that  whole 
rock  systems  can  become  partially  open  (Cumming  and  Scott 
1976,  Hoffman  1971,  Brooks  1980,  Olszewski  1980)  or  even 
completely  reset  during  metamorphic  events  (Taylor  1975, 
Weber  et  al.  1975). 

Although  it  has  been  demonstrated  that  common 
strontium  can  also  become  mobile  under  metamorphic 
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conditions  (Baadsgaard  and  van  Breemen  1970),  it  seems 
reasonable  to  assume  that  under  most  conditions  its 
mobility  on  a  whole  rock  scale  would  be  negligible 
compared  to  rubidium  and  radiogenic  components  of 
strontium.  The  bulk  of  the  normal  strontium  component 
is  usually  distributed  within  the  feldspars  where  it 
is  tightly  held  in  the  host  lattice,  even  at  elevated 
temperatures.  The  radiogenic  strontium  resulting  from 
decay  of  8 ;Rb  is  usually  present  in  a  mineral  lattice 
which  is  a  relatively  good  host  for  rubidium  but  not 
for  radiogenic  strontium.  Radiogenic  strontium,  being 
poorly  diadochic  in  this  lattice,  will  have  a  much 
greater  tendency  to  diffuse. 

It  is  clear  that  if  the  scale  of  strontium  isotope 
re-equilibration  and  rubidium  migration  is  much  less 
than  the  size  of  the  collected  and  analysed  specimens, 
the  whole  rock  Rb-Sr  data  should  fall  on  an  isochron 
recording  the  initial  age  of  the  sample.  According  to 
Roddick  and  Compston  (1977) ,  the  formation  of  a  large 
rock  unit  may  be  imagined  as  consisting  of  a  set  or  sub¬ 
systems,  each  of  which  is  large  enough  to  have  an  Rb-Sr 
ratio  equivalent  to  that  of  the  rock  unit  as  a  whole. 

If  the  range  of  migration  is  larger  than  the  scale  of 
one  such  subsystem,  then  the  whole  rock  Rb-Sr  data 
should  fall  on  an  isochron  recording  the  age  of  the 
metamorphic  event.  If,  on  the  other  hand,  the  range  of 
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migration  is  less  than  one  such  subsystem  but  larger 
than  that  of  the  analysed  sample,  scatter  will  be 
observed  in  the  standard  isochron  diagram. 

It  was  said  at  the  beginning  of  this  chapter  that 
the  size  of  the  collected  specimens  was  on  order  of 
1-2  lb.,  except  the  samples  collected  from  the  central 
part  of  the  Pegasus  Lake  anticline,  which  were  about 
twice  as  large.  It  might  not  be  a  coincidence  that  data 
obtained  on  these  samples  show  a  tendency  to  align 
themselves  along  the  trend  line,  suggesting  the  age 
of  initial  emplacement  of  the  Charlebois  Lake  complex 
precursors  (figure  14).  The  scatter  of  these  data 
suggests  that  the  range  of  migration  was  probably  not 
much  larger  than  the  size  of  the  collected  specimens. 

The  other  smaller  samples  collected  from  the  outer  areas 
of  the  anticline,  those  of  granodior itic  granofels  and 
those  from  the  "Old  Camp  Area",  exhibit  scatter  and  some 
of  them  show  a  tendency  to  align  themselves  along  the 
isochrons  recording  the  age  of  the  metamorphic  event. 
This  behaviour  suggests  two  alternatives.  First,  the 
range  of  migration  and  re-equilibration  was  about  the 
same  for  the  centrally  located  samples  from  the  Pegasus 
Lake  anticline  and  for  the  samples  from  the  other  areas. 
However,  the  samples  from  the  other  areas  were  smaller, 
resulting  in  the  metamorphic  age  being  "more  visible". 
The  other  alternative  is  that  although  the  samples  were 
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smaller,  the  range  of  migration  in  the  "Old  Camp  Area" 
and  other  areas  close  to  and  including  the  granodiori tic 
granofels  was  much  larger  than  in  the  core  of  the  Pegasus 
Lake  anticline.  If  this  were  the  case,  uranium  mineraliza¬ 
tion,  rubidium  enrichment  in  the  granodioritic  granofels 
and  formation  of  migmatites  could  have  at  least  partially 
been  caused  by  this  local,  small-scale  remobilization  and 
migration  of  the  material  during  the  Hudsonian  meta- 
morphic  event,  as  suggested  by  Morra  (1977).  Morra ' s 
suggestion  (1977)  that  sedimentary  precursors  of  the 
granodioritic  granofels,  migmatites  and  granitic  gneiss 
from  the  outer  margins  of  the  anticline  have  been 
originally  enriched  in  uranium  and  potassium  (therefore 
rubidium  also)  appears  to  be  supported  by  the  lead-lead 
and  Rb.-Sr  data  presented  so  far.  It  is  clear  from  lead- 
lead  and  Rb-Sr  data  that  lead  and  strontium  of  the  grano¬ 
dioritic  granofels  have  evolved  in  high  y  and  high  Rb-Sr 
environments  during  premetamorphic  time,  while  the  same 
elements  of  the  granitic  gneisses  from  the  core  of  the 
Pegasus  Lake  anticline  have  evolved  in  lower  y  and  Rb-Sr 
environments . 

Rb-Sr  data  obtained  on  six  whole  rock  samples  of  the 
5  m  wide  late  felsic  dyke  which  crosscuts  other  rock  units 
in  the  "Old  Camp  Area"  are  presented  in  Table  VI  and  Figure 
17.  The  data  exhibit  a  range  of  e;Rb/86Sr  ratios  from 
around  6.0  to  over  30.0.  All  of  the  points  scatter  rather 
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Figure  17.  Rb-Sr  diagram  for  the  whole  rock  samples  of  the  pink  felsic  dyke  displays  general  scatter  of 
data  points  around  the  1300  m.y.  reference  isochron  and  high  (87Sr/86Sr)  ratio.  For  identi¬ 
fication  of  the  specimens  refer  to  Figure  4  and  Table  VI .  Data  ellipses  are  drawn  at  4a 
error  level. 


badly  so  that  no  reasonable  isochron  may  be  computed. 
An  approximate  line  drawn  through  the  scattered  points 
gives  an  approximate  age  of  about  1300  m.y.  and  an 
initial  87Sr/86Sr  ratio  of  about  0.730.  Although  this 
age  does  not  have  a  statistical  certainty  because  of 
very  wide  scatter,  it  is  strikingly  similar  to  the  age 
obtained  from  lead-lead  data  on  the  same  samples. 
Contamination  by  the  metasedimentary  rocks  (Baadsgaard 
and  Godfrey  1972)  and  remelting  of  the  rocks  from  the 
lower  crust  might  be  responsible  for  the  scatter  and 
high  apparent  initial  87Sr/86Sr  ratio. 
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CHAPTER  IV 


SUMMARY  AND  CONCLUSIONS 


The  study  was  undertaken  to  establish  the  age  and  the 
history  of  the  area  and  to  possibly  glean  some  information 
regarding  the  relationship  between  the  individual  members 
of  the  complex  and  the  origin  of  the  U-Pb  system  of  the 
uranium  mineralized  rock  units.  It  became  apparent,  early 
in  this  study,  that  interpretation  of  the  geochronological 
data  for  this  area  would  be  complicated  by  the  open  system 
behaviour  and  limited  migration  of  elements  during  the  last 
metamorphic  event.  For  that  reason  a  complete  set  of  lead- 
lead  and  rubidium-strontium  whole  rock  analyses  was  per¬ 
formed  on  the  same  specimens.  Radiogenic  minerals  such  as 
uraninite  and  mineral  fractions  rich  in  monazite  from  the 
granodior itic  granofels  and  two  migmatite  samples  (i.e. 
samples  with  the  highest  206Pb/204Pb  ratios)  were  separated 
and  analysed  for  lead,  uranium  and  thorium  isotopes.  Six 
feldspar  separates  from  the  least  radioactive  samples, 
that  is, samples  with  the  lowest  206Pb/204Pb  ratios,  were 
analysed  for  lead  isotopes. 

The  geochronological  data  obtained  during  this  study 
was  correlated  with  the  available  geological  evidence  and 
summarized  in  Table  VII. 

The  uranium-lead  and  lead-lead  results  seem  to  have 
been  more  useful  than  the  rubidium-strontium  results, 
which  exhibited  general  scatter  and  provided  only  slight 
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suggestions  of  the  important  events  in  the  history  of 
these  rocks.  The  uranium-lead  data  for  uraninites  and 
monazite-rich  fractions  show  clear  indications  of  a 
relatively  slow  and  gradual  weakening  of  the  Hudsonian 
metamorphic  event  beginning,  in  this  area,  about  1797  m.y. 
ago  with  the  closing  of  the  U-Pb  system  of  monazites  and 
possibly  other  high-grade  minerals.  U-Pb  data  obtained 
on  the  uraninite  separates  suggest  that  the  end  of  the 
Hudsonian  orogeny  in  the  Charlebois  Lake  area  was  marked 
by  the  closing  of  the  U-Pb  system  of  uraninite  1750  m.y. 
ago.  The  ages  of  368  m.y.  and  55  m.y.  obtained  from  the 
apparent  lower  concordia  intercepts  of  the  uraninite  and 
monazite  discordia  lines,  respectively,  are  believed  to 
be  spurious  since  discordance  is  the  result  of  continuous 
lead  diffusion  from  these  minerals. 

The  overall  history  of  the  Charlebois  Lake  complex 
is  most  completely  illustrated  by  the  lead-lead  and,  to 
some  extent,  supported  by  the  rubidium-strontium  results. 
Whole  rock  and  feldspar  lead-lead  data  indicated  that  lead 
in  this  complex  approximates  the  evolution  according  to 
the  three  stage  evolution  model  with  f^  =  const.  (Gale 
and  Mussett  1973) .  It  is  suggested  that  the  original 
emplacement  (i.e.  formation  of  the  sedimentary  precursors) 
of  the  Charlebois  Lake  complex  could  have  taken  place  as 
early  as  2470  m.y.  ago  from  an  environment  that  had 
y  =  9.85.  This  value  is  slightly  higher  than  Stacey  and 
Kramers  (1975)  estimate  of  the  earth's  y2  value.  The 
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rubidium-strontium  data  obtained  on  the  specimens  from 
the  core  of  the  Pegasus  Lake  anticline  show  vague  signs 
that  0.701  might  have  been  the  initial  average  87Sr/86Sr 
ratio  of  the  sedimentary  precursors  some  2400  m.y.  ago. 

From  the  overall  average  87Rb/86Sr  of  0.412  for  the 
Charlebois  Lake  complex,  the  mean  87Sr/86Sr  ratio  of 
around  0.705  at  the  end  of  the  Hudsonian  metamorphic 
event  and  assuming  that  the  initial  87Sr/86Sr  could  not 
have  been  very  much  different  from  0.701,  it  was  estimated 
that  the  age  of  the  sedimentary  precursors  must  be  around 
2440  m.y.  It  therefore  was  concluded  from  all  these 
indications  that  the  Charlebois  Lake  complex  was  probably 
formed  2400-2470  m.y.  ago. 

Besides  the  lead-uranium  data  of  the  radiogenic 
minerals,  the  Hudsonian  metamorphic  event  was  well  defined 
by  the  parallel  lead  distribution  of  the  lead-lead  data 
with  slope  ages  ranging  from  1741  m.y.  to  1775  m.y.  The 
rubidium-strontium  data  also  provided  some  indications  of 
this  event,  expressed  in  trend  lines  formed  by  grano- 
dioritic  granofels  samples,  specimens  from  the  outer 
margins  of  the  Pegasus  Lake  anticline  and  the  lower  dis¬ 
tribution  boundary  of  the  "Old  Camp  Area"  data. 

Morra  (1977)  observes  that  there  is  no  evidence  of 
large  scale  migration  of  material  in  the  area,  and  suggests 
that  uranium  mineralization  is  probably  a  product  of 
locally  remobilized  uranium  (and  possibly  thorium)  from 
the  uranium  and  thorium-rich  sedimentary  layers  of  the 
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arkosic  arenites.  These  suggestions  are  supported  by 
the  lead-lead  and  rubidium-strontium  data  presented  in 
this  study.  The  distribution  of  the  lead-lead  and  to 
some  extent  the  rubidium-strontium  data,  indicated  that 
direct  precursors  of  the  mineralized  granodior itic 
granofels  had  evolved  in  high  y,  i.e.  a  high  uranium  and 
possibly  high  rubidium  environment  during  premetamorphic 
time.  Both  lead-lead  and  rubidium-strontium  data  of  the 
mineralized  granodior it ic  granofels  represent  the  more 
radiogenic  limits  of  the  distribution  of  the  geochrono- 
logical  data.  The  other  samples  exhibit  a  gradual 
decrease  in  radiogenic  component  from  mineralized  grano- 
dioritic  granofels,  which  is  adjacent  to  the  calc-silicate 
layer  towards  the  more  massive  granitic  gneiss  at  the 
base  of  the  metasedimentary  succession.  One  may  speculate 
that  the  increase  in  pressure  during  compaction  of  loose 
sediments  "squeezed  out"  fluids  which  served  to  facilitate 
migration  of  uranium  upwards.  The  precursors  of  the  calc- 
silicates  (calcareous  and  tuffaceous  sediments)  acted  as 
an  impermeable  barrier  to  the  migrating  solution,  result¬ 
ing  in  the  relatively  higher  concentration  of  uranium 
adjacent  to  this  layer.  Later,  under  metamorphic  conditions 
uranium  was  locally  remobilized  producing  the  observed 
radioactive  anomalies  and  the  distribution  of  lead-lead 
data  presented  in  Figure  10. 

Due  to  the  insufficient  number  of  collected  and 
analysed  samples,  the  history  of  the  late  felsic  dykes 


and  the  irregular  pegmatitic  bodies  has  not  been  well 
determined.  Lead-lead  and  rubidium-strontium  data  indicate 
that  this  rock  unit  was  emplaced  in  the  pos t-Hudsonian 
time,  possibly  some  time  between  1500  m.y.  and  1100  m.y. 
ago.  An  extensive  geochronological  investigation  using 
lead-lead  and  rubidium-strontium  methods  of  this  and  other 
similar  late  rock  units  from  the  surrounding  areas  and 
larger  Churchill  province  is  strongly  suggested.  The 
benefit  of  such  a  study  would  be  a  better  understanding 
of  the  history  and  events  taking  place  in  the  Churchill 
province  during  this  time  period. 

This  study,  like  many  others  in  the  literature,  has 
demonstrated  the  necessity  of  the  multimethod  approach  in 
geochronological  investigations  of  disturbed  rock  systems. 
Although  it  was  shown  that  the  whole  rock  Pb-Pb  system 
preserves  the  history  record  well,  and  possibly  better 
than  the  Rb-Sr  system,  the  use  of  two  or  more  methods 
generally  reveals  more  details  and  allows  more  confidence 
in  making  geochrological  conclusions  and  assumptions. 
Whenever  enough  material  is  available,  the  use  of  combined 
Pb-Pb,  Rb-Sr  and  U-Pb  dating  technique  on  whole  rock  and 
mineral  separates  is  strongly  suggested. 
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APPENDIX  I 


THEORY  OF  U-Th-Pb,  Pb~Pb  AND  Rb~Sr  SYSTEMATICS 
U-Th-Pb  Systemat ics 

Three  basic  radioactive  decay  systems  are  used  in 
U-Th-Pb  and  Pb-Pb  age  dating  techniques 

2  3  8u^2  0  6pb  +  8  i(  He  +  6S~ 

2  3  5  u^2  0  7pb  +  7  4He  +  48“ 

a 

2  3  2Th^2  0  8pb  +  64He  +  4  8~  # 

The  relationship  of  the  molar  concentration  of  the 
daughter  isotope  (corrected  for  initial  daughter  contamina¬ 
tion)  to  the  molar  concentration  of  the  parent  isotope, 
measured  in  the  sample  at  present  (t=0)  is  given  by  the 
following  general  equation: 

D  =  P [ exp (At)  -  1] 

where  D  is  the  molar  concentration  of  the  daughter 
isotope,  P  is  the  molar  concentration  of  the  parent 
isotope,  presently  measured,  A  is  radioactive  decay 
constant  and  t  is  the  time  when  the  isotopic  system 
became  closed  to  migration  of  either  parent  or  daughter 
isotope  across  the  boundaries  of  the  system. 

Radioactive  decay  constants,  adopted  throughout  this 
study,  are  (as  suggested  by  Steiger  and  Jager  1977): 
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Decay 

constant 

of 

2  3  8  U  : 

A  =  0. 155.125x10  yr 

Decay 

constant 

of 

2  3  5U  : 

A*  =  0 . 98485xl0~9yr~1 

Decay 

constant 

of 

2  3  2  Th . 

A  * '  =  0. 049475xl0“9yr“1 

It  is  immediately  obvious  that  knowing  the  amounts 
of  parent  and  daughter  isotopes  at  present,  the  time 
when  the  system  became  closed  can  be  calculated  from  the 
relationship : 

t  =  1/A  An  +  1]  . 

Applying  this  equation  to  the  three  decay  systems 
mentioned  at  the  beginning,  three  independent  ages  can 
be  obtained. 

Since  the  present  relative  isotopic  proportions  of 

natural  uranium  are  known  ((2j5U/2  l8U)  ,  =  137.88) 

present 

and  are  assumed  to  have  varied  over  geologic  time  due  to 
the  radioactive  decay  only,  it  is  possible  to  obtain  a 
fourth  useful  date,  i.e.  707Pb/2GbPb  age. 

By  a  comparison  of  these  ages  (i.e.  measured 
daughter  to  parent  ratios  of  the  sample  system  being 
investigated)  it  is  possible  to  derive  not  only  the  age 
of  the  system  but  some  information  regarding  the  history 
of  the  isotopes  and  therefore  the  sample  system  as  well. 
If  all  the  calculated  ages  are  ,the  same  within  analytical 
error,  the  phase  being  dated  has  remained  closed  to 
changes  in  the  parent-daughter  system  and  such  dates  are 


termed  "concordant  dates". 


Theoretical  U-Pb  systematics  are  graphically 
expressed  in  Figure  18.  It  shows  that  the  mineral  or 
other  investigated  phase  which  has  not  been  subjected  to 
the  lead  loss  or  uranium  gain  (or  lead  gain  or  uranium 
loss)  will  have  concordant  207Pb/23"u,  206Pb/“38U  and 
2  0  7  Pb//  2  0  6  Pb  ages  and  that  data  will  be  represented  by  a 
point  that  lies  on  the  curved  line  called  "concordia" . 

If  the  dated  phase  has  been  subjected  to  the  lead 
loss  (or  uranium  gain)  during  the  short  period  of  time, 
recently  or  any  time  in  the  past  after  the  original 
closure  of  the  system,  obtained  ages  will  not  agree, 
within  the  set  error  limits,  i.e.  they  will  be  "discor¬ 
dant".  U-Pb  data  obtained  on  a  suite  of  samples  that 
has  experienced  such  episodic  daughter  losses  (or  parent 
gains)  whether  recently  or  in  a  distant  past,  will 
theoretically  give  a  chord  of  points  that  determines  the 
straight  line  which  intersects  the  concordia  curve  at 
two  places.  The  upper  intersection  represents  the  time 
of  original  formation  of  the  phase  and  the  lower  inter¬ 
section,  the  time  of  the  episodic  event  which  caused  the 
disturbance  in  the  system. 

If  it  is  assumed  that  the  lead  continuously  diffuse 
from  the  system  over  the  entire  history  of  the  phase, 
measured  Pb-U  ratios  on  a  set  of  samples  with  this  type 
of  history  will  fall  along  the  continuous  diffusion  line 
in  the  concordia  plot  (Figure  18)  . 
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Figure  18.  Concordia  diagram  showing  various  types  of  discordance  all  starting  at  2700 
m.y.  I:  episodic  lead  loss  line.  II:  continuous  lead  diffusion  line. 

Ill:  recent  lead  loss  (from  concordant  mineral)  line. 
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A  combination  of  these  basic  types  of  parent- 
daughter  histories  can  greatly  complicate  the  interpre¬ 
tation  of  the  obtained  data.  A  more  complete  discussion 
of  the  U-Th-Pb  systematics,  along  with  the  set  of  impor¬ 
tant  references,  is  given  in  Faure  (1977). 

Pb-Pb  Systematics 

Detailed  treatment  of  Pb-Pb  systematics  is  covered 
elsewhere  in  the  literature  (Faure  1977,  Gale  and  Mussett 
1973,  Kanasewich  1968,  and  others)  so  that  only  a  brief 
overview  of  the  theory  as  applied  to  the  whole  rock  system 
is  given  here. 

Ratios  of  the  radioactive  isotopes  2UGPb,  2 0 7Pb  and 
208Pb  (generated  according  to  the  radioactive  decay 
schemes,  mentioned  in  the  previous  section)  to  the  non- 
radiogenic  2  0  4  Pb  isotope  are  utilized  for  interpretation. 

A  more  useful  way  of  presenting  the  data  is  on  a  207Pb/ 

2  0  4Pb  vs.  206Pb/204Pb  rather  than  on  a  208Pb/204Pb  vs. 
206Pb/204Pb  diagram.  The  reason  for  this  is  that  uranium 
isotopes  (whose  products  are  207Pb  and  206Pb)  cannot  be 
significantly  separated  during  the  normal  geochemical 
reaction,  whereas  uranium  and  thorium  (whose  decay  product 
is  2  0  8  Pb)  may  be . 

The  interpretation  of  the  ]read  isotope  distribution 
in  a  suite  of  rock  samples  basically  consists  of  fitting 
it  to  the  model  which  would  be  described  by  a  set  of 
episodic  changes  in  U~Pb  ratio  at  some  discrete  time 
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intervals  in  the  past. 

The  simplest  case  of  the  lead  isotope  evolution  would 
be  if  the  lead  evolved  by  a  simple  radioactive  decay  of 
uranium  to  lead  over  geologic  time.  This  evolution  is 
described  by  the  following  equations: 

X]  =  ao  +  ~  exp  (  Xt  1 )  ] 

yx  =  b0  +  (y/137.88)  [exp(A’to)  -  expU't^j 

where  x^  =  206Pb/204Pb  and  y^  =  207Pb/204Pb,  both  at  the 

time  t_ 7  a  =  206Pb/204Pb  and  b  =  207Pb/204Pb,  both  at 
1  o  o 

the  time  t  ;  p  =  238U/204Pb  at  present  time  (t=0) ,  A  and 
A'  are  radioactive  decay  constants  with  the  values  given 
in  the  previous  section.  The  time  t  is  measured  back  into 
time  from  the  present.  When  plotted  on  the  x-y  (207Pb/ 
2°4pb  versus  206Pb/204Pb)  diagram  these  two  equations  are 
represented  by  a  curve  which  is  termed  "lead  evolution 
growth  curve"  (Figure  19).  The  values  of  a^  and  bQ  are 
lead  isotopic  ratios  in  the  triollite  phase  of  the  Canyon 
Diablo  meteorite  i.e.,  aQ=9.307,  bQ=10.294,  as  suggested  by 
Tatsumoto  et  al.  (1973).  It  is  assumed  that  isotopic 
composition  of  the  Pb  in  the  earth  was  equal  to  a^  and  bQ 
at  a  time  which  is  referred  to  as  the  age  of  the  earth. 

This  term,  which  is  not  an  easily  measurable  value,  can 
be  estimated  in  two  ways.  First,  based  on  the  slope  of 
the  line  which  is  determined  by  the  207Pb/204Pb  and 
206Pb/204Pb  ratios  of  various  meteorites,  and  second 


1 


101 


L 


O 

* 

10 

H 


o 

o 

CO 

tH 

tH 

'1— 1 

q  d 

f  02 


/ 

/ 

/ 


4GZ 


o 

o> 


o 

co 


Figure  19.  Pb  growth  curve  diagram.  y  =  9.85,  a  =  9.346,  b  -  10.218  and  T  —  4570  m.y. 
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based  on  the  isotopic  ratios  of  conformable  lead  ore 
deposits  which  appear  to  fall  on  a  single  stage  growth 
curve.  Their  ages,  as  determined  by  their  position  on 
the  lead  growth  curve,  appear  to  be  somewhat  younger 
than  their  geologic  ages  determined  by  other  methods. 
Discrepancies  become  more  pronounced  with  younger  ages, 
suggesting  that  the  single  stage  lead  evolution  growth 
curve  fits  the  history  of  these  ores  only  as  a  first 
approximation.  More  sophisticated  models  of  mantle 
crustal  evolution  that  take  into  account  possible  cata¬ 
strophic  changes  (Stacey  and  Kramers  1975)  or  continuous 
changes  (dimming  and  Richards  1975)  in  the  U-Pb  system 
have  been  devised.  From  a  model  with  a  continuous  change 
in  y  Cumming  and  Richards  estimated  the  age  of  the  earth 
to  be  4509  m.y.  while  Stacey  and  Kramers  adopted  (from 
Tatsumoto  et  al.  1973)  4570  m.y.  for  the  age  of  the 

earth  and  also  suggested  global  scale  fractionation  of 
the  U-Pb  system  at  3700  m.y.  ago. 

In  a  closed  system  the  only  changes  taking  place  in 
the  U-Pb  ratio  are  discrete  and  episodic.  It  may  be  said 
that  the  suite  of  "U-Pb  subsystems"  which  was  separated 
from  the  global  y^  system  at  time  t^ ,  with  Pb-Pb  isotopic 
composition  x  ^  ^  and  y^  and  then  each  subsystem  evolved 
in  separate  y0^  environments  until  time  t ^  (when  all  the 
uranium  was  removed)  will  distribute  themselves  along  the 
straight  line  of  slope: 
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y(2)  “  y(1)  l  expiX't^  -  exp(A't2) 

x ( 2 )  ~  ~  137.88  exp(At^)  -  exp(At2) 

If  the  system  evolved  until  the  present,  then  t2=0, 
and  the  above  slope  equation  may  be  solved  directly  for 
t^  (Fig.  20).  This  is  a  two  stage  model.  The  line  passes 
through  the  primary  growth  curve  of  the  y^  environment  at 
t=0  and  t=t^.  Any  of  these  subsystems  can  be  characterized 
by  two  parameters:  and  f^  =  y^/y.^.  The  parameter  f 

is  the  fractionation  coefficient  and,  providing  there  is 
no  change  in  the  isotopic  composition  of  uranium  or  lead 
at  the  time  of  differentiation  (t^ ) ,  it  is  a  direct  measure 
of  the  U-Pb  ratio  immediately  before  and  after  the  diffe¬ 
rentiation  . 

The  general  three  stage  model  history  involves 
evolution  in  the  constant  y  environment  from  the  formation 
of  the  earth  to  the  time  ( t^ )  of  the  first  global  fractiona¬ 
tion  of  the  system  into  subsystems  of  different  ^(i) 
environments.  Each  second  stage  subsystem  evolves  inde¬ 
pendently  until  the  new  U-Pb  fractionation  of  these 
individual  subsystems  takes  place  at  time  into  the  third 
stage  environments  which  are  characterized  by  the  P1 3  (  -l  j  ) 
values.  From  then  until  the  present,  (t=0)  each  third 
stage  subsystem  remains  closed  ahd  evolution  continues  by 
the  decay  of  uranium  into  lead  only. 

Each  episodic  event  is  characterized  by  a  correspond¬ 
ing  fractionation  coefficient:  at  t^  by  f^  =  y1/y9;  and 
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Figure  20.  Pb  isotopic  evolution  diagram  for  the  samples  which  underwent  two-stage 
evolution  history.  Initial  y  =  9.85,  Tq  =  4570  m.y.,  t^  =  4000  m.y., 

a  =  9.346  and  b  =  10.218. 


The  resulting  general  distribution  of  the  points  in 
the  2  0  7Pb/ 2  0  4  Pb  vs.  206Pb/204Pb  diagram  will  be  governed 
by  the  degree  of  fractionation  at  t^  and  t?.  The  most 
common  case  will  be  when  f^  and  are  both  variable  and 
the  points  distribute  themselves  in  a  fan-shaped  fashion 
with  the  apex  denoting  the  lead  composition  which  the 
system  possessed  at  t^ .  The  lower  boundary  of  the  dis¬ 
tribution  of  data  points  will  have  a*  slope  age  approaching 
t^  while  the  upper  boundary  will,  in  an  ideal  situation, 
intersect  the  primary  growth  curve  at  t-j  and  t0.  This 
behaviour  can  be  most  easily  understood  by  considering 
Figure  21.  Lead  evolution  according  to  the  three 
stage  model  generally  results  in  scattered  points  on 
an  X-Y  diagram.  Only  in  the  special  cases  of  f^  = 
constant  or  f^  =  constant  are  linear  plots  obtained. 

The  f  =  const. /f 2  =  variable  model  is  a  generally  more 
plausible  model  from  a  geological  standpoint  and  theore¬ 
tically  results  in  a  family  of  parallel  straight  lines 
with  a  constant  slope  which  depends  only  on  t^,  (time  of 
the  second  fractionation)  and  not  on  t^  (time  of 
the  first  U-Pb  fractionation) .  If  f,  =  1  there  is  no 
fractionation  of  U-Pb  system  from  y  at  time  t^,  therefore 
the  evolution  history  reduces  to  a  two-stage  model.  If 
f  becomes  very  large  (i.e.  large  uranium  loss  or  lead 
gain  (with  the  same  isotopic  composition)  occurs  at  time 
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t^) ,  the  limiting  straight  line  will  intersect  the  primary 
growth  curve  at  t^.  Low  values  of  imply  uranium  gain 
or  lead  loss  (or  both)  occurred  at  the  first  fractionation 
event  ( t ^ )  .  In  the  case  of  f1  =  variable/f„  =  constant 
model,  theoretically  a  family  of  straight  lines  all  passing 
through  the  same  point,  i.e.  t^  of  y^  growth  curve  will 
result.  The  slopes  and  upper  intercepts  with  the  primary 
growth  curve  do  not  have  any  real  meaning.  This  model  is 
less  geologically  plausible  since  it,  requires  constant 
fractional  change  for  variable  y  environments  during  the 
second  fractionation  event  at  time  t 9 .  It  is  also  some¬ 
what  contradictory  in  the  case  of  the  grand  system  remaining 
closed.  The  model  allows  for  lead  loss  at  the  second 
fractionation  event  but  not  for  a  gain  of  lead  with 
different  isotopic  compositions.  If  the  global  system 
remains  closed  and  some  of  the  subsystems  with  different 
y's  lose  some  of  their  lead  (which  has  a  different  iso¬ 
topic  composition) ,  other  subsystems  will  be  liable  to 
gain  that  lead,  violating  the  conditions  for  obtaining  a 
rectilinear  plot  for  a  suite  of  samples  analysed  at  present. 
A  model  consisting  of  more  than  three  episodic  stages  would 
produce  an  overall  distribution  similar  to  the  three  stage 
model . 
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of  the  molar  concentration  of  the  daughter  isotope  to 
the  molar  concentration  of  the  parent  isotope,  measured 
in  the  sample  at  present  (t=0),  is  given  by  the  following 
equation : 

8  7Sr  =  8  7Rb [ exp ( At )  -  1]  +  8 7SrQ 

where  0  ;Rb  is  the  molar  concentration  of  the  8  7Rb  isotope 
and  8  ;Sr  is  the  molar  concentration  of  the  8 'Sr  isotope, 
measured  at  present.  A  is  the  radioactive  decay  constant 
of  8 7Rb  =  1.42x10  "^yr  ^  (as  suggested  by  Steiger  and 
Jager  (1977)),  t  is  the  time  of  the  initial  isotopic 
homogenization  measured  positively  into  the  past  and  87SrQ 
is  the  amount  of  8 ;Sr  originally  present. 

Because  of  the  presence  of  the  other  Sr  isotopes 
and  the  convenience  of  handling  isotopic  ratios,  the  above 
equation  is  modified  to  include  concentration  of  the  86Sr 
isotope  into: 

8  7Sr/ 8  6  Sr  =  8 7Pb/86Sr [exp(At)  -  1]  +  ( 8 7Sr/ 8 6 Sr )  . 

This  parametric  equation  (of  the  type  y  =  mx  +  n) 
contains  two  different  types  of  information  that  could  be 
determined  from  two  or  more  cogenetic  samples,  i.e.  from 
the  slope  and  initial  (8;Sr/86Sr)  .ratio  from  the  inter¬ 
cept.  The  model  age  determination  from  the  single  rock 
or  mineral  specimen  requires  an  assumption  about  the 
initial  strontium  isotopic  composition  ( 8 7Sr/ 8 6 Sr )  . 
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The  data  from  a  suite  of  samples  which  have  been 
isotopically  homogenized  at  the  same  time  (t)  in  the  past 
may  be  evaluated  in  two  ways.  The  more  popular  way  is  by 
plotting  87Sr/86Sr  versus  8/Rb/86Sr  values.  A  oogenetic 
suite  of  samples  which  have  been  completely  homogenized 
at  the  time  the  system  was  last  disturbed  should  form  a 
straight  line  known  as  an  isochron  (Figure  22) . 

Another  way  of  representing  the  data  is  by  following 
the  evolution  of  87Sr/36Sr  with  timejf  i.e.  by  plotting 
87Sr/86Sr  versus  time  (Figure  23).  This  is  achieved  by 
calculating  the  value  of  87Sr/86Sr  at  a  time  approximately 
equal  to  the  time  of  formation  and  joining  this  point  to 
the  presently  measured  8/Sr/86Sr  ratio.  By  doing  this 
for  the  oogenetic  suite  of  samples,  an  array  of  lines 
will  be  obtained.  The  slopes  of  the  lines  are  propor¬ 
tional  to  the  8/Rb/86Sr  ratio  of  each  sample  since  its 
formation.  All  the  lines  intersect  at  the  time  equal  to 
the  time  of  isotopic  homogenization  and  at  the  87Sr/86Sr 
ratio  that  the  suite  had  at  that  time. 

The  only  possible  advantage  of  this  treatment  over 
the  standard  isochron  diagram  is  that  it  occasionally 
gives  a  clearer  view  of  the  evolution  of  the  systems 
which  have  been  disturbed  after  the  initial  homogeniza¬ 
tion.  All  other  information  obtained  from  this  diagram 
is  also  obtainable  from  the  isochron  diagram. 

The  samples  selected  for  dating  may  be  a  set  of 
whole  rocks  selected  from  the  single  formation  to  give  a 
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Figure  22.  Rb-Sr  diagram  for  the  set  of  samples  which  were  emplaced  1800  m.y.  ggo 
with  an  initial  87Sr/86Sr  =  0.703. 


Ill 


Figure  23.  Sr  growth  lines  diagram  for  the  1800  m.y.  old  suite  of  oogenetic  samples. 


broad  range  of  Rb-Sr  values,  or  they  could  be  mineral 
separates  from  a  single  rock  such  as  low-to-medium  Rb-Sr 
K-feldspars  or  apatites  and  high  Rb-Sr  micas. 

The  usefulness  of  the  Rb-Sr  dating  method  is  entirely 
dependent  on  the  complete  homogenization  of  the  Rb  and  Sr 
isotopes  at  some  time  in  the  past  and  upon  minerals  and 
rock  remaining  closed  systems  to  the  migration  of  Rb  and 
Sr  from  that  time  to  the  present. 

If  the  samples  to  be  dated  are  /crystallized  from 
magma  and  there  was  no  subsequent  disturbance  of  the 
system,  through  thermal  or  other  events,  the  mineral 
separates  and  whole  rock  samples  should  give  the  same 
true  age  of  the  rock. 

If  the  samples  have  been  subjected  to  metamorphism 
of  sufficient  intensity  to  recrystallize  some  of  the 
minerals  or  make  them  open  to  diffusion  of  Rb  and  Sr, 

Rb-Sr  data  of  minerals  will  record  the  time  of  metamor¬ 
phism  and  the  Rb-Sr  data  of  whole  rocks  the  time  of 
initial  homogenization.  This  will  hold  true  providing 
the  range  of  Rb  and/or  Sr  migration  has  been  limited  to 
an  area  smaller  than  the  size  of  the  whole  rock. 

A  complete  review  of  Rb-Sr  systematics  is  given  in 
Faure  and  Powell,  1972  and  Faure,  1977. 
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APPENDIX  II 


ANALYTICAL  METHODS 

General  Laboratory  Procedures 

During  the  course  of  this  work  most  of  the  separa¬ 
tions,  filament  loadings  and  other  operations  were  carried 
out  in  the  filtered  air  laboratory  or  in  a  laminar  flow 
hood  equipped  with  an  HEPA  100+  filter.  The  water  used 
has  been  purified  by  passing  through  a  large  mixed  bed  ion 

t 

exchange  column  and  then  distilling  it  twice  in  the  sub- 
boiling  stills.  The  first  distillation  was  done  in  a 
glass  still;  water  was  fed  into  the  smaller  all-quartz 
still  where  a  second  distillation  was  completed.  Nitric 
and  hydrochloric  acid,  of  reagent  grade,  were  purified, 
first  by  distillation  in  an  all-quartz  sub-boiling  still 
and  then  in  a  slow  "two  teflon  bottle"  vapour  distillation 
still  (Mattinson  1972).  Hydrofluoric  acid  and  hydrobromic 
acid  of  ARI STAR  grade  (BDH  Chemicals  Ltd.)  were  further 
purified  by  distillation  in  "two  teflon  bottle"  stills. 

The  sulphuric  acid  and  perchloric  acid  used  were  of 
ULTREX  grade  supplied  by  J.T.  Baker  Ltd. 

All  of  the  isotope  tracers  used  in  this  study  were 
obtained  from  the  National  Bureau  of  Standards  or  Oak 
Ridge  Laboratories  and  further  calibrated  and  prepared  by 
Dr.  H.  Baadsgaard  of  this  department. 

Lead-Lead  Analyses 

One  to  two  grams  of  finely-ground  whole  rock  or 
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K-feldspar  sample  was  weighed  out  in  a  150  ml  teflon 
beaker.  Ten  ml  of  concentrated  HF  and  5  ml  of  concen¬ 
trated  HNO  were  added  to  the  beaker  and  stirred.  The 
beaker  was  then  covered  with  a  teflon  "watch  glass"  and 
left  overnight  at  approximately  80°C.  The  cover  was 
removed,  the  solution  evaporated  to  dryness  and  the 
residue  "baked"  for  one  hour  at  105°C  to  remove  excess 
HF .  The  residue  was  moistened  with  a  few  drops  of  con¬ 
centrated  HNO ^  and  "baked"  again  for  one-half  hour  to 
convert  fluoride  salts  to  nitrate  salts.  This  step  was 
repeated  three  times. 

Ten  mis  of  8N  HNO^  was  added  to  the  residue,  the 
beaker  was  covered  and  heated  at  80 °C  for  two  hours.  The 
solution  and  gelatenous  precipitate  was  transferred  to  a 
15  ml  centrifuge  tube  and  the  gel  was  centrifuged  off. 

The  solution  was  poured  back  into  the  teflon  beaker  and 
evaporated  to  dryness. 

A  few  drops  of  2N  HC£  and  5  ml  of  IN  HBr  were  added 
to  the  dried  residue  (for  larger  samples  more  HBr  is 
used),  the  beaker  was  covered  and  left  overnight.  The 
residual  gel  (if  any)  was  centrifuged  off  before  anion 
exchange  separation. 

The  quartz  columns  used  for  anion  exchange  separa¬ 
tion  of  lead  were  filled  with  CGA  540-X8,  100-200  mesh, 
anion  exchange  resin,  supplied  by  J.T.  Baker  Ltd.  Later 
on  during  this  study,  the  resin  was  replaced  with  AG  1-X8, 
100-200  mesh,  anion  exchange  resin,  supplied  by  BIO-RAD 
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Laboratories . 

A  resin  bed  of  0.4  cm  x  12.5  cm  was  prepared  by 
passing  30  ml  of  ^0,  30  ml  of  9N  HC£  and  30  ml  of  H^O, 
and  then  equilibrated  with  30  ml  of  IN  HBr .  The  sample 
solution  was  passed  through  the  column,  followed  by  15  ml 
of  IN  HBr  and  3  ml  of  2N  HC£.  The  lead  was  eluted  with 
10  ml  of  9N  HC £  and  collected  in  the  small  teflon  beaker. 
While  the  lead  solution  was  being  evaporated  to  dryness, 
the  column  was  prepared  for  HC£  ion  .exchange  purification 
of  lead  by  passing  through  it  an  additional  20  ml  of 
9N  HC£  followed  by  30  ml  of  H^O .  The  resin  bed  was  set 
with  20  ml  of  2N  HC£. 

The  residue  containing  lead  in  the  small  teflon 
beaker  was  dissolved  in  0.5  ml  of  2N  HC£  and  carefully, 
without  disturbing  the  resin  bed,  loaded  on  the  column. 
The  solution  was  washed  in  with  3x2  ml  of  2N  HC£  and  the 
lead  was  eluted  in  the  fraction  between  6  and  12  ml  of 
2N  HC  £ . 

This  fraction  was  evaporated  to  dryness  and  lead 
chloride  was  redissolved  in  a  few  drops  of  Ho0  (preferred 
concentration  at  this  point  is  approximately  500  yg 
Pb/1  ml  of  solution) . 

Lead  was  loaded  on  a  single,  zone  refined,  rhenium 
filament,  using  the  silica  gel  technique:  first  the  drop 
of  the  aqueous  lead  solution  was  gently  evaporated  onto 
the  rhenium  filament  (the  usual  amount  of  lead  on  the 
filament  is  2-3  ygm) .  A  large  drop  of  silica  gel  slurry 


* 

in  0 . 8N  H^PO^  was  loaded  on  top  and  the  reaction  observed 
through  the  microscope.  The  current  passing  through  the 
filament  was  increased  to  1.2  amps.  The  excess  water 
evaporated,  the  cloudy  gel  solution  began  to  liquify  and 
clear  up.  At  1.4  amps  the  drop  on  the  filament  was 
completely  clear  and  at  1.6  amps  lead-phosphosilicate 
started  precipitating.  At  around  2.0  amps  phosphoric 
acid  was  fumed  off.  The  filament  was  brought,  briefly, 
to  "dull  red"  by  increasing  the  current  to  approximately 
2.5  amps. 

Maximum  lead  blank,  for  up  to  2.5  gm  of  rock  sample 
was  measured  to  be  not  more  than  15  nanograms,  which  is 
generally  well  below  0.5  per  mil  of  total  lead  in  the 
sample  and  is  considered  to  have  a  negligible  effect  on 
the  isotopic  ratio  in  this  study. 

A  filament  with  lead  sample  was  introduced  into 
ALDERMASTON  MICROMASS  MM30  solid  source  mass  spectro¬ 
meter.  After  the  pump-down,  the  temperature  of  the 
filament  was  raised  to  1400°C  (as  determined  by  the 
optical  pyrometer) .  The  initially  erratic  lead  beam 
stabilizes  usually  after  30  minutes.  The  mass  spectro¬ 
meter  was  operated  in  the  peak  switching  mode  and  with 
the  Faraday  cup  collection  system.  The  acceleration 
voltage  used  was  8  kv  and  the  data  collection  and 
reduction  were  performed  on  line  by  an  HP9825  computer. 
The  measuring  error  of  the  instrument  was  better  than 
0.15  per  mil  (standard  error  of  the  weighted  mean  ratio) 
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J  0  4Pb  measuring  error:  0.31°/oo 

Overall  la  errors  assigned  to  this  type  of  ratio  in  this  study  are: 
2  0  6  pb/  2  0  4pb  =  o . 3  6  °/0  o  2  0  7Pb/  2  0  4Pb  =  0 . 50°/oo  208Pb/204Pb  =  0.70 
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The  fractionation  correction  factors  for  measured  lead 
.ratios,  the  overall  error  assigned  to  the  studied  samples 
and  the  reproducibility  of  measurements  were  determined 
by  repetitive  analyses  of  NBS  SRM-981  common  lead  standard. 
The  obtained  results  are  given  in  Table  VIII. 

Lead- Uranium-Thorium  Analyses 

The  samples  analysed  for  lead-uranium  or  lead- 
ur an ium- thorium  were  uraninites  and  monazite-r ich  mineral 
fractions . 

The  minerals  were  originally  separated  from  their 
respective  powdered  whole  rock  samples  utilizing  magnetic 
and  heavy  liquid  mineral  separation.  Three  of  the  samples 
contained  uraninite  grains  large  enough  to  be  picked  under 
the  microscope.  The  fourth  sample  contained  some  uraninite 
intergrown  with  monazite  and  the  fifth  sample  did  not  have 
any  visible  uraninite.  Intergrown  uraninite  was  leached 
with  cold  4H  HNO ^  for  three  hours. 

From  the  remainder  of  the  heavy  minerals,  impurities 
such  as  pryite,  molybdenite,  occasional  biotite  flakes 
and  other  mineral  grains  were  removed  under  the  micro¬ 
scope  so  that  a  mineral  concentrate  containing  more  than 
50%  monazite  was  obtained.  The  total  amount  of  the 
samples  available  ranged  from  less  than  0.1  mg  for  tiny 
uraninite  cubes  to  up  to  30  mg  for  large  monazite  con¬ 
centrates  . 

The  uraninite  samples  were  dissolved  in  4N  HNO^  and 
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monazite-rich  fractions  (after  washing  in  warm  2N  HNO^) 
and  were  dissolved  in  5  ml  of  concentrated  HF  and  2  ml 
of  concentrated  HNO ^  in  a  teflon  decomposition  vessel 
similar  to  that  of  Krogh  (1973).  The  monazite  solution 
was  evaporated  to  dryness  and  the  residue  converted  to 
nitrate  form  by  treating  it  with  a  few  drops  of  concen¬ 
trated  HNO^  and  then  it  was  redissolved  in  4N  HNO^. 

Both  uraninite  and  monazite  solutions  were  divided  into 
two  aliquots.  The  larger  ones  were  .used  for  Pb-IR 
determination  and  the  smaller  ones  for  Pb-U-Th  concen¬ 
tration  determination.  To  this  aliquot  208Pb,  ^ 3 5U  and 

230Th  tracers  of  known  amounts  were  added  and  the  solution 
evaporated  to  dryness  and  redissolved  in  a  drop  of  2N  HC£ 
and  2  ml  of  IN  HBr. 

Lead  from  both  Pb-IR  and  Pb-U-Th-ID  aliquots  was 
purified  in  the  same  manner  as  described  previously  in 
the  section  on  whole  rock  lead  separations.  The  mixture 
of  Pb,  U  and  Th  from  the  isotope  dilution  aliquot  was 
loaded  on  the  HBr  column,  the  column  was  eluted  with  IN 
HBr  and  2N  HC£  (lead  was  "stripped"  later  with  9N  HC£) 
and  the  collected  solution  containing  U  and  Th  was 
evaporated  to  dryness.  The  residue  was  carefully  converted 
to  nitrate  form  with  a  few  drops  of  concentrated  HNO^  and 
again  evaporated  to  dryness. 

In  the  meantime,  a  glass  column  was  filled  with 
AG1-X8, 100-200  mesh,  anion  exchange  resin  in  nitrate  form. 

A  resin  bed  of  0.8  cm  x  12.5  cm  was  prepared  by  passing 
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50  ml  of  H ,?0 ,  30  ml  of  7N  HNO^  and  50  ml  of  H^O  and 
finally  equilibrated  with  30  ml  of  7N  HNO^* 

The  dry  residue  containing  uranium  and  thorium  was 
picked  up  in  a  few  drops  of  7N  HNO^  and  carefully  loaded 
on  the  resin  bed.  The  column  was  eluted  with  6><2  ml  of 
7N  HNO^,  after  which  uranium  and  thorium  were  "stripped" 
off  the  column  with  20  ml  of  This  solution  was 

evaporated  to  dryness,  the  uranium  and  thorium  were 
redissolved  in  a  drop  of  H „0  and  the  drop  was  carefully 
evaporated  on  the  rhenium  side  filament  of  the  bead 
equipped  with  a  double  filament.  On  the  center  filament 
the  corresponding  Pb~ID  was  loaded,  so  that  all  three 
concentration  measurements  were  performed  using  one  double 
filament  bead. 

The  bead  was  introduced  in  the  source  chamber  of  the 
MM3Q  mass  spectrometer.  After  the  lead  analysis  (from 
the  center  filament)  was  completed,  the  center  filament 
temperature  was  raised  to  2130°C.  The  current  flowing 
through  the  side  filament  was  increased  until  a  steady 
uranium  beam  was  obtained.  Uranium  peak  switching  data 
collection  and  reduction  were  done  by  an  HP9825  computer. 

In  order  to  obtain  a  thorium  beam,  uranium  was 
"burned  off"  by  increasing  the  side  filament  current  by 
about  one  additional  ampere.  The  ionization  of  thorium 
began  when  the  center  filament  temperature  was  raised  to 
2230°C.  Data  collection  starts  as  soon  as  the  thorium 
beam  reaches  20%  intensity  on  a  "GAINX100"  scale. 
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Overall  errors  on  207Pb/235U,  206Pb/238U  and 

208Pb/232Th  are  relatively  difficult  to  assess,  due  to 
the  many  possible  sources  of  random  and  systematic  error, 
such  as  mass  discrimination,  blank,  etc. 

Radiogenic  lead  measurements  have  not  been  corrected 
for  mass  discrimination  since  it  is  negligible  for  highly 
radiogenic  samples.  The  presently  assumed  level  of 
accuracy  for  the  Pb/U  and  Pb/Th  ratios  in  the  laboratory 
is  better  than  1%  at  two  sigma.  These  error  limits  have 
been  assigned  to  all  the  measurements  in  this  study  and 
therefore  have  been  used  for  line  fittings  and  age  cal¬ 
culations  . 

Rubidium-Strontium  Analyses 

Powdered  whole  rock  samples  were  analysed  by  X-ray 
fluorescence  spectrometry- to  give  an  estimate  of  Rb  and 
Sr  concentrations.  From  these  results  enough  sample  was 
used  to  provide  about  15-20  yg  of  strontium  and  15-20  yg 
of  rubidium  for  analyses. 

Rubidium  separation: 

Enough  sample  was  weighed  out  into  a  teflon  beaker 
to  provide  approximately  15  yg  of  rubidium,  8 7Rb  tracer 
was  added  along  with  10.0  ml  of  concentrated  HF  and  a 
few  drops  of  concentrated  HN00.  The  beaker  was  covered 
with  teflon  "watch  glass"  and  heated  for  a  few  hours  at 
approximately  80°C.  After  the  decomposition  of  the 
sample  was  completed,  the  solution  was  transferred  to 
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the  platinum  dish,  5  drops  of  concentrated  H  SO  were 
added  and  the  solution  was  evaporated  to  dryness.  The 
temperature  was  raised,  excess  H^SO^  was  fumed  off  and 
the  sample  was  ignited  at  800°C  for  30  minutes  with  the 
Bunsen  burner.  When  cooled,  the  residue  was  leached  with 
a  few  drops  of  Ho0  and  the  solution  (0.25-0.5  ml  in  total) 
was  transferred  to  a  3  ml  centrifuge  tube.  Several  drops 
of  pure  concentrated  HCfO^  were  added,  resulting  in  the 
precipitate  of  HC£0^  and  coprecipitated  RbC&O  .  The  tube 
was  centrifuged  and  the  supernate  was  poured  off.  The 
walls  of  the  tube  were  then  rinsed  carefully  with  water. 

A  drop  or  two  of  water  was  added  to  the  precipitate  so 
that  some  of  it  dissolves.  A  small  amount  of  this  solu¬ 
tion  was  carefully  evaporated,  and  then  fused  on  the 
rhenium  side  filament. 

Strontium  separation: 

The  powdered  rock  sample  was  weighed  out  into  a 
150  ml  teflon  beaker,  the  sides  of  the  beaker  were  rinsed 
with  several  ml  of  water  to  bring  down  material  which 
might  have  clung  electrostatically  to  the  walls.  A  pre¬ 
weighed  84Sr  tracer  was  quantitatively  added  to  the 
sample  in  the  teflon  beaker.  The  beaker  was  then  left 
uncovered  at  low  heat  so  that  the  water  would  slowly 
evaporate.  When  the  residue  was  dry,  10  ml  of  concen¬ 
trated  HF  and  5  ml  of  concentrated  HNG^  were  added  to  the 
beaker  and  after  stirring  gently  the  beaker  was  covered 


with  teflon  "watch  glass"  and  left  overnight  at  approxi¬ 
mately  80 °C.  The  cover  was  then  removed,  the  solution 
evaporated  to  dryness  and  the  residue  was  "baked"  at 
105°C  for  one  hour.  The  fluoride  salts  were  converted 
to  nitrate  form  by  moistening  the  residue  with  concen¬ 
trated  HNO^  and  evaporating  it  to  dryness.  This  step  was 
repeated  three  times.  Ten  mis  of  8N  UNO  were  added  to 
the  beaker  and  heated  gently  for  a  few  hours,  stirring 
occasionally  to  redissolve  the  residue.  The  solution 
was  evaporated  to  half  of  its  original  volume,  transferred 
to  the  centrifuge  tube,  the  undissolved  silica  gel  (if 
any)  was  centrifuged  off.  The  supernate  was  temporarily 
transferred  to  its  teflon  beaker,  the  gel  from  the  tube 
was  discarded,  the  tube  was  rinsed  and  the  supernate  was 
returned  to  the  centrifuge  tube. 

Three  drops  of  concentrated  Ba(NO^)„j  solution  were 
added  to  the  sample  solution,  stirred,  and  the  total  volume 
made  up  to  10  ml  with  concentrated  HNO^.  The  sides  of  the 
centrifuge  tube  were  rubbed  with  the  teflon  rod  to  initiate 
precipitation  of  fine  BaCNO^)^  crystals,  which  coprecipitat 
strontium  ions.  After  standing  for  two  hours,  the  Ba{NO^)^ 
was  centrifuged  off  and  the  supernate  was  discarded. 

During  the  latter  part  of  this  study  Rb-Sr  analyses 
were  performed  using  mixed  Rb-Sr  "spiking"  techniques. 

This  involves  adding  mixed  tracer  to  the  powdered  rock 
sample  and  following  the  strontium  separation  procedure  to 
the  point  of  centrifuging  off  Ba(NO^)^. 


The  supernate 
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containing  rubidium  is  transferred  to  a  platinum  dish, 

5  drops  of  H  SO.  are  added  and  the  rubidium  separation 
Z  4 

procedure  is  followed. 

The  Ba(NO^)^  precipitate  was  washed  with  concentrated 
HNO_  and  then  recrystallized  by  dissolving  it  first  in  a 
minimum  amount  of  H^O  and  reprecipitating  it  by  adding 
8  ml  of  concentrated  HNO^  and  rubbing  the  walls  of  the 
tube  with  a  teflon  rod.  After  standing  for  two  hours, 
the  precipitate  was  centrifuged  off,'  the  supernate  was 
discarded,  the  tube  was  washed  carefully  with  2N  HC£  so 
that  the  precipitate  is  not  disturbed,  and  finally  the 
precipitate  was  dissolved  in  approximately  0.5  ml  of 
2N  HC£. 

A  quartz  cation  exchange  column  was  prepared  to 
separate  barium  and  strontium.  An  0.4  cm  x  16  cm 
column  bed  (approx.  2  ml  by  volume)  of  DOWEX  50W-X12, 
100-200  mesh  cation  exchange  resin  was  prepared  by 
passing  30  ml  of  H^O,  30  ml  of  9N  HC£,  30  ml  of  H^O  and 
finally  30  ml  of  2N  HC£.  The  solution  containing  barium 
and  strontium  was  carefully  loaded  on  the  column  without 
disturbing  the  resin  bed.  The  column  was  eluted  with 
2N  HC£  and  strontium  was  collected  in  the  eluted  fraction 
between  12  ml  and  20  ml. 

The  strontium  solution  was  evaporated  to  dryness, 
SrC£2  was  picked  up  in  a  small  drop  of  water  and  gently 
evaporated  on  a  rhenium  side  filament  of  a  double 


filament . 
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Mass  spectrometr ic  analyses  were  performed  using  an 
ALDERMASTON  MM30  mass  spectrometer.  Later,  rubidium 
measurements  were  made  using  a  home-built  12  inch,  90° 
magnet  sector,  solid  source  mass  spectrometer  with  an  on 
line  TI  980A  minicomputer.  This  instrument  was  operated 
in  the  peak  switching  mode,  with  5  KV  acceleration  voltage 
and  a  Faraday  cup  collection  system. 

Strontium  and  rubidium  blanks  for  the  initial  samples 
which  were  done  by  spiking  separate  ,aliquots  of  samples, 
were  measured  to  be  a  2  ng  strontium  blank  and  a  4  ng 
rubidium  blank.  Blank  measurements  on  the  mixed  spike 
procedure  produced  a  5  ng  strontium  blank  and  an  8  ng 
rubidium  blank.  In  both  cases  the  effect  of  the  blank 
was  considered  negligible  and  no  blank  corrections  were 
applied  to  the  results. 

87Sr/86Sr  measuring  precision  and  overall  errors, 
assigned  to  the  8/Sr/06Sr  ratios  of  the  analysed  samples, 
were  determined  by  repetitive  analyses  of  NBS  SRM-987 
strontium  carbonate  standard  and  by  comparing  the  obtained 
results  with  the  NBS  suggested  values.  It  was  found  that 
the  results  obtained  in  this  study  were  slightly  higher 
than  values  suggested  by  NBS  (Table  IX) .  Since  the 
discrepancies  are  very  small,  no  corrections  of  the 
measured  87Sr/86Sr  ratios  were  made. 

Rubidium  fractionation  correction  coefficient  and 


reproducibility  were  determined  by  repetitive  measure¬ 
ments  of  isotopic  ratios  performed  on  the  NBS  SRM-984 


1 


rubidium  chloride  standard  and  on  rubidium  extracted 
from  the  common  granite  sample  (Table  IX) . 

The  individual  errors  for  fl7Rb/06Sr  are  somewhat 
difficult  to  determine  precisely  because  of  the  many 
possible  sources  of  systematic  error  such  as  mass  dis¬ 
crimination,  blank,  etc. 

The  accepted  laboratory  accuracy  for  e7Rb/86Sr 


during  this  study  was  1%  at  two  sigma  level. 


TABLE  IX 


Isotopic  Analyses  of  Rubidium 
and  Strontium  Standards 


NBS  suggested  value: 
8  7Rb/ 8  5 Rb  =  0 .3856 

Measured  8  7  Rb/ 8  b  Rb 


NBS  suggested  value: 
8  7Sr/ 8  6  Sr  =  0.71014 

.0002 

Measured  8 7 Sr/ 8  6  S r 


0.38513  SRM-984 

0.38533 

0.38488 

0.38500 

0.38429 

0.38523  Corrrnon  Granite  Sample 

0.38401 

0.38433 

0.38477 


0.71034  SRM-987 

0.71019 
'0.71033 
0.71024 
0.71032 
0.71011 
0.71010 
0.71007 
0.71007 


Mean:  0.38477  0.71021 

±.00048  ±.00011 

la  error:  1.25°/00  la  error:  0.15°/oo 

Rubidium  fractionation  correction  factor: 


( 8  ?Rb/ 8  5Rb) measured  _  0 . 3 8477 
( 8 yRb/ H  bRb)  true  0.3  8  56 


0.99772 


V 


